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ABSTRACT 

Correlations between intrinsic shear and the density field on large scales, a poten- 
tially important contaminant for cosmic shear surveys, have been robustly detected 
at low redshifts with bright galaxies in SDSS data. Here we present a more detailed 
characterization of this effect, which can cause anti-correlations between gravitational 
lensing shear and intrinsic ellipticity (GI correlations). This measurement uses 36 278 
Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) spectro- 
scopic sample with 0.15 < z < 0.35, split by redshift and luminosity; 7 758 LRGs from 
the 2dF-SDSS LRG and QSO (2SLAQ) Survey at 0.4 < z < 0.8; and a variety of 
other SDSS samples from previous, related work. We find > 3cr detections of the effect 
on large scales (up to 60 hr 1 Mpc) for all galaxy subsamples within the SDSS LRG 
sample; for the 2SLAQ sample, we find a 2a detection for a bright subsample, and 
no detection for a fainter subsample. Fitting formulae are provided for the scaling of 
the GI correlations with luminosity, transverse separation, and redshift (for which the 
2SLAQ sample, while small, provides crucial constraints due to its longer baseline in 
redshift). We estimate contamination in the measurement of as for future cosmic shear 
surveys on the basis of the fitted dependence of GI correlations on galaxy properties. 
We find contamination to the power spectrum ranging from —1.5 (optimistic) to —33 
per cent (pessimistic) for a toy cosmic shear survey using all galaxies to a depth of 
B = 24 using scales I « 500, though the central value of predicted contamination is 
—6.5 per cent. This corresponds to a bias in tx 8 of Aa s = —0.004 (optimistic), —0.02 
(central), or —0.10 (pessimistic). We provide a prescription for inclusion of this error 
in cosmological parameter estimation codes. The principal uncertainty is in the treat- 
ment of the L ^ L* blue galaxies, for which we have no detection of the GI signal, but 
which could dominate the GI contamination if their GI amplitude is near our upper 
limits. Characterization of the tidal alignments of these galaxies, especially at redshifts 
relevant for cosmic shear, should be a high priority for the cosmic shear community. 

Key words: cosmology: observations - gravitational lensing - large-scale structure 
of Universe. 
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Weak gravitational lensing has emerged in recent years as 
a powerful tool for probing the distribution of matter in 
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the universe (see e.g. the review bv lRefregien l2003h. Lens- 
ing directly traces the matter distribution, and thus is less 
sensitive to modeling of baryonic physics than other cos- 
mological probes. Because of the relative cleanness of the 
theory, several groups are engaged in a programme to mea- 
sure cosmological parameters using the two-point function 
of weak lensing she ar. The field has grown rapidly , with 
the first detections jBacon et aL 2000l; Kaiser et al.l |2000| ; 



IVan Waerbeke et aLlfcood ; IWittman et al.ll2000h giving way 
to more robust results with smaller statistic al errors and 
better control of obse r vational systematics dBrown et al 



20031; iJarvis et al1l2003l ; iHevmans et~ai1l2005l ; iMassev et al 



20051 ) . Several of the more recent results provid e competitive 
cosmological constraints dHoekstra et al ] |2006l ; IJarvis et alj 
l200d : ISembolini et alj|2006h . 

In principle, because the cosmic shear signal has the 
power to constrain the matter power spectrum amplitude 
at a particular redshift, it may be a powerful tool to con- 
strain models for dark energy. Splitting source samples into 
redshift slices, and auto- and cross-correlating each slice 
against itself and against the other slices can constrain the 
growth of perturbations as a function of time. Indeed, lens- 
ing tomography was shown t o put sign i ficant constra i nts on 
the dark energy paramet ers l|Hull2002l ; iHutererlliool ; llshakl 
120051 ; ISong fc Knoxll2005l ) and to have th e potential to test 
gravitational phy sics on very large scales (jlshak et al.ll2006l ; 
iKnox et alj|200q ). Therefore, many future surveys are being 
planned to have sufficiently high statistical power to mea- 
sure dark energy using cosmic shear. 

The cosmological weak lensing signal is small, with typ- 
ical shear of order 10 -2 or less, and some of the more am- 
bitious projects under consideration will require measure- 
ment of the signal to very high fractional accuracy (< 1 
per cent). This means that even very small systematic er- 
rors can have a significant influence on cosmic shear anal- 
yses. One of these systematics is the problem of measur- 
ing shape in the presence of an instrumental point spread 
function (PSF), which can easily introduce spurious ellip- 
ticity correlations between different galaxies and dilute ex- 
isting correlations by making the galaxies appear rounder. 



There is now a substantial literature on the sub 


ect of 


PSF-induced systematics llKaiser 2000; 


Erben et al. 


2001; 


Bernstein & Jarvis 2002: Hirata & Seliak 


2003a; Ishak et al. 


2004; Vale et al.l 20041; IVan Waerbeke et al.l 2005T). includ- 



ing a d etailed comparison of different measureme nt algo- 
rithms l|Hevmans et al.l l2006al ; IMassev et all 120071 ). A sec- 
ond problem is the determination of the redshift distribu- 
tion of source galaxies, AN/dz, which enters into the equa- 
tion for the cosmic shear power spectrum. There are sev- 
eral ways in which this can be addresse d, such as direct 
spectroscopic measurement of some fields ([Bernstein fc Jainl 
120041 ; llshak fc Hiratall2005l: IMandelbaum et al.ll2005h . self- 
calibration ( Huterer et al. T2006t ) , and possib ly in the future 
by la rge-scale surveys in the H I 21 cm line (Rawlings et al. 
l2004h . Yet a third potential difficulty is the uncertainty in 
the theory: the cosmic shear programme places very de- 
manding requirements on the accuracy of iV-body simula- 
tions, and the effec t s of baryons will be important on sm all 
scales l| White! l2004l ; |Zhan fe Knoxll2004 ITing et al.ll20o!l ). 

Most of the above-described systematics are technical 
in nature, and in principle can be overcome or at least sup- 
pressed with improved modeling and data reduction tech- 



niques. Unfortunately, there is another possible systematic 
in weak lensing, namely intrinsic alignments of the source 
galaxy population. It has been recognized for some time 
that if the intrinsic ellipticities of different galaxies are cor- 
related, then this can result in a spurious contribution to the 
"shear" powe r spectrum. Several th e oretical model s were 



constructed JCroft fc Metzlei 2000l; iHeavens et al 



ICatelan et~aill200ll ; ICrittenden et alj l200ll : |jindl2002l ) and 



2000; 



the intrinsic alignment correlation function was measured 
by several authors using lo w-redshift surveys (|Brown et al.l 
120021 ; IHevmans et alj|2004h . The severity of this problem is 
lessened by the fact that such intrinsic ellipticity-intrinsic 
ellipticity (II) correlations can only exist if the two source 
galaxies are near each other in physical space, whereas the 
shape correlations induced by cosmic shear exist between 
any two galaxies that are near each other in angular coordi- 
nates, even if they have a large line-of-sight separation. Thus 
several authors proposed that galaxy pairs at neighbouring 
redshi fts be down-weighted or remove d in cosmic shear anal- 
yses (iKing fc Schneiderl |2002i , 120031 ; IHevmans fc Heavens! 
120031 ; iTakada fc White! 120041 '). However, it was later found 
that cosmic shear is subject to another type of intrinsic 
alignment contamination, in which the intrinsic ellipticity of 
a nearby galaxy correlates with the quadrupole of the den- 
sity field surrounding it. Since this density field produces a 
lensing shear on more distant galaxies, it is possible to have 
a gravitational lensing-intrinsic ellipticity (GI) correlation. 
The GI correlation can in principle exist between galaxies 
at very different redshifts, and therefore it cannot be elim- 
inated by considering o nly shear cross-correla tions between 
different redshift slices ijHirata fc Seliakll2004r ). 

The cosmic shear community has proposed several 
methods for addressing the GI pr oblem, such as marginaliza- 
t ion over parameterized models |Kindl2005l ; iBridle fc Kind 
l2007h and removal based on the redshift dependence of the 
signal (|Hirata fc Seliak 2004). At the same time, it would 
be useful to know how large a GI signal to expect, including 
dependence on the redshift and on the source galaxy lumi- 
nosity, colour, environment, etc. Models of GI correlation, 
either empirical or theoretical, would provide an indication 
of how much of the contaminant needs to be removed in or- 
der to realize the full potential of cosmic shear, and would 
help guide removal strategies (e.g. rejecting certain types of 
galaxies with the worst GI contamination) . At present, there 
are s everal analytical dHui fc Zhandl2002l; iHirata fc Seliakl 
l2004h and simulation-based l|Hevmans et al.ll2006bl ) models, 
while observationally the GI correlatio n function relevant 
for lens ing has been measured only by IMandelbaum et al.l 
(|2006ah . There is a much larger literature on galaxy align- 
ments that uses statistics other than the GI correlation func- 
tion, or_^neaj^es_rorrelafe 

pie. IPen et al.ll20qol; ILee fc Penll200ll; lBernstein fc Norben 



20021 ; ILee fc Pen! |2002|; IHirata et al.l 120041; iNavarro et al 
20041; lAgustsson fc Brainerdl 120061: lAzzaro et all 120071 : 



Donoso. O'Mill. fc Lamba s 2006; Yang et al.l I2006T ). These 



studies are useful for constraining theoretical models, but 
cannot readily be turned into quantitative empirical esti- 
mates of contamination to cosmic shear on large scales. 

In a previous paper (Mandel baum et al.ll2006ah we ob- 
tained a detection of the large-scale density-ellipticity corre- 
lation and presented a preliminary analysis of the expected 
contamination to cosmic shear results. This paper extends 
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that original work by providing a phenomenological charac- 
terization of the intrinsic alignment 2-point correlation func- 
tions that is as comprehensive as can be obtained with exist- 
ing data. We use a variety of galaxy samples from the SDSS 
Main spectroscopic sample, the SDSS spectroscopic Lumi- 
nous Red Galaxy (LRG) sample, and the 2dF-SDSS LRG 
and QSO (2SLAQ) Survey. In particular, we explore the lu- 
minosity, colour, environment, and redshift dependence of 
the signal. This is an important step in constructing and 
validating models (either empirical or simulation-based) of 
intrinsic alignments. In this paper we will obtain scaling re- 
lations and construct an empirical model for the intrinsic 
alignment 2-point functions for red galaxies; comparison to 
simulations and consideration of 3-point functions will be 
addressed in a future work. 

We begin in Section[5]by reviewing briefly the formalism 
for describing the cosmic shear and intrinsic alignment cor- 
relation functions. The sources of data used for this work are 
described in Section [3] The methodology used for the mea- 
surement is detailed in Section [4] The results obtained and 
systematics tests are discussed in Section [S] The bias of the 
LRG samples, needed to estimate the density, is estimated in 
Section [6] We present power- law fits for the GI correlations 
in Section [7] and compare against theoretical predictions in 
Section [8] Using these fitting formulae, we estimate pro- 
jected contamination for measurements of as with current 
and future cosmic shear surveys in Section [§] we discuss how 
cosmic shear analyses can marginalize over our estimates 
of GI contamination. We conclude in Section 1101 Two ap- 
pendices are included: Appendix [X] provides our correlation 
function data, and Appen dixlBlpr ovid es the calculations for 



terms, 
C EE (a(3) 



C 



EE,GG 
I 



(af3) + Cf EJI (a0) 



+Cf E ' GI {al3). (1) 

In iMandelbaum et al.l (|2006a) we presented the Limber in- 
tegrals that allow us to determine each of these quantities 
in terms of the matter power spectrum and intrinsic align- 
ments power spectrum. Of greatest interest to us is the GI 
contamination term, 



C 



f («/3) = J ^f a (r)W (r)P ml Ik-- 
0), 



1 + 1/2 



+ (a - 



(2) 



where m is the comoving distance to the horizon, f a (r) is 
the comoving distance distribution of the galaxies in sample 
a, and 



W a (r) 



f a (r')dr'. 



(3) 



(The generaliz ation of these equations to curved universes 
can be found in lMandelbaum et alj|2006al .) The power spec- 
trum P s -i(k) is defined as follows. If one chooses any two 
points in the SDSS survey, their separation in redshift space 
can then be identified by the transverse separation r p and 
the radial redshift space separation FlQ The + and x com- 
ponents of the shear are measured with respect to the axis 
connecting the two galaxies (i.e. positive + shear is radial, 
whereas negative + shear is tangential) . Then one can write 
the density-intrinsic shear correlation in Fourier space as 



function data, and Ap pendixlBlprovides the calculations tor r 

converting the lHevmans et alT(|2006bh simulation results to P s^' = ~ 2yr J &+( r J" n ) M kr p) r v dr p dn , 



(4) 



the variables used in this paper. 



2 FORMALISM 



The form alism for the analysis o f the lensing shear two-point 
function (|Miralda-Escudelll99ll ) and the intrinsic alignment 
contamination is well developed. We will briefly summarize 
the important equations here, and then define some new 
variables that relate to observables in g alaxy surveys. Our 
not ation is consistent with tha t of lllirata fc Seliak (|2004T ) 
and IMandelbaum et all l|2006ah . 



where £s+(r p ,Tl) is the correlation function between 5 and 
7+. Here y 1 is the intrinsic shear weighted by the density of 
galaxies, i.e. 7 1 = (1+Sg)"/ 1 , and 5 and 5 g refer to the matter 
and galaxy overdensities, respectively. It is often convenient 
to do the projection through the radial direction, 



(•») 



w s+ (r p ) = / £ s +(r p , II)dII. 



A similar set of equations can be written for the intrinsic- 
intrinsic terms; se e e.g. Eqs. (5), (7), (9), and (10) of 
IMandelbaum et all (|2006al ). 



2.1 Correlation functions 

The observed shear 7 of a galaxy is a sum of two compo- 
nents: the gravitational lensing-induced shear 7°, and the 
"intrinsic shear" 7 7 , which includes any non- lensing shear, 
typically due to local tidal fields. Therefore the _E-mode 
shear power spectrum between any two redshift bins a and 
(3 is the sum of the gravitational lensing power spectrum 
(GG), the intrinsic-intrinsic, and the gravitational-intrinsic 



1 Technically the "environment dependence of the 2-point func- 
tion" is some combination of higher-order statistics, not a 2-point 
function. However it can be treated by the same methods used to 
understand colour or luminosity dependence, so we investigate it 
here. 



2.2 Cosmology dependence 

Here we note the cosmological model and units adopted 
for this work and the effect of changing it. Pair separa- 
tions are measured in comoving h^ 1 Mpc (where Ho = 
lOOfo kms -1 Mpc -1 ), with the angular diameter distance 
computed in a spatially flat ACDM cosmology with Q m — 
0.3. The correlation function measurement depends on Q m 
via the determination of both r p and n. If fi m ^= 0.3, then 
our calculated r p and n are in error via the equations 



r p (calc) _ Da(z;Q,„ 



r p (true) Da (z; 0.3) 



and 



n(calc) _ H(z; 0.3) 
n(truc) ~ H(z;Q m ) 



(6) 



2 The redshift space separation is frequently denoted tt; we use II 
to avoid confusion since the number it = 3.14... appears frequently 
in this paper. 
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Note that since angular diameter distance Da is mea- 
sured in h^ 1 Mpc and Hubble rate H(z) is measured in 
kms -1 (h^ 1 Mpc)^ 1 , there is no dependence of these equa- 
tions on h. The errors on r v and II are simple proportional- 
ities and hence when we do power-law fits to the projected 
correlation function, w g +(r p ) = Ar p , there is no effect on 
the power-law index a. The amplitude will however be in 
error by 



(7) 



A(calc) _ H(z;0.3) l D A (z-n m ) 
A(true) ~~ H(z;Q m ) |_ D^z; 0.3) 

This equation should be used when interpreting our results 
in other cosmologies. It should be noted however that the ef- 
fect is rather small: for a = —0.8 (roughly what we observe) 
and z = 0.6 (the cosmology matters most at high redshift), 
the ratio J 4(calc)/ J 4(true) = 1.05 for Q. m = 0.2 and 0.96 for 
Q, m = 0.4. Since this encompas ses the range of fi m values 
from recent determinations (e.g. lSpergel et al.ll2007l ). the er- 
ror in our results due to uncertainty in Q m is thus at the 
~ 5 per cent level, which is negligible compared to our final 
errors. Therefore we have not included it explicitly in our 
error analysis, and have not attempted to re-measure the 
correlation function for different cosmologies. 



3 DATA 

Much of th e data used here are obtained from the SDSS. 
The SDSS (|York et al.l l2000h is an ongoing survey to im- 
age ~ tv steradians of the sky, and fo llow up ~ 10 6 



of the detected objects spectroscopically ( Eisenstein et al] 
l200ll ; IStrauss et alj|2002l ; iRichards etafl 120021 ) ■ The imag- 



ing is carried out by drift-scanning the sky in ph oto- 



Ivezic et al. 2004). in 



metric conditions dHogg et all l200ll; . 
five bands (ugriz) (|Fukugita et al.lll996l ; ISmith et alj|2002h 
a specially designed wide-field camera (|Gunn et al.l 



1998). These imaging data are the source of the LSS sam- 
ple that we use in this paper. In additio n, objects are 
target ed for spectroscopy using these data (|Blanton et al.l 
2003a:) and are observe d with a 640-fibr e spectrograph 
on the same telescope l|Gunn et al] [2006). All of these 
data are processed by completely automated pipelines 
that detect and measure photometric pro perties of ob- 
jects , and astrom e trical ly calibrate the data l|Lupton et al] 
120011 ; iPier et al l 120031 ; iTucker et all I2006T). The SPSS 



n^haTsjTr^ ioT^t a releases kought on et al 



Abazaiian et al.l I200I |2004 [20051 : iFinkbeiner et al 



2002 



2004 



Adelman-McCarthv et all 2006 a. b). This analysis uses the 



spectroscopically observed galaxies in the Value-Added 
Gala xy Catalog, LSS sample 14 (VAGC; iBlanton et al] 
120051 ). comprising 3423 square degrees with SDSS spectro- 
scopic coverage, as well as photometric catalogs to be de- 
scribed below covering a larger area. 

To determine density-shear correlation functions, one 
needs a sample of galaxies that traces the intrinsic shear, 
and another that traces the density field. Unlike in 
iMandelbaum et alj l|2006al ). for this paper we use a vari- 
ety of galaxy samples spanning redshifts z — 0.05 — 0.8. The 
sources of these samples are the SDSS Main spectroscopic 
sample (z = 0.05 — 0.2), the SDSS spectroscopic LRG sample 
[z — 0.16 — 0.36), and the 2SLAQ spectroscopic LRG sample 
[z — 0.4 — 0.8). In the following subsections we describe each 



of these samples, and Table \T\ summarizes their characteris- 
tics. In each of these three redshift regions, the full sample 
is used to trace the density field, but we break the galax- 
ies down into subsamples based on luminosity and (for the 
Main sample) colour to trace the intrinsic shear field. This 
enables us to explore the possibility that different types of 
galaxies show different intrinsic alignment signals. 

3.1 SDSS Main subsamples 

The first samples of galax ies we discuss are those used in 
IMandelbaum et al.1 i|2006al ). subsamples of the SDSS Main 
spectroscopic sample divided by luminosity and other prop- 
erties. The GI correlation models in this paper are split into 
early and late types, so we have presented here the GI signals 
for the Main galaxies split both by luminosity, and by colour 
as "blue" or "red." (The differences bet ween the Main sam- 
ple an alysis of this section and that of iMandelbaum et al] 
l2006al are the inclusion of the colour separator, and the use of 
the full Main sample rather than the luminosity subsamples 
to trace the density field.) For this work, we use the galax- 
ies in L3 through L6, one magnitude fainter than L* (L3) 
through two magnitudes brighter (L6). The sample proper- 
ties were described in full in that paper; for this work, we 
mention only that the luminosities described are Petrosian 
magnitudes, extinctio n corrected using reddening maps from 
ISchlegel et al.l (|l998Tl with the extinc tion-to-reddening ra- 
tios given in IStoughton et al.l (|2002T l. and fc-corrected to 
z = 0.1 using KCOR RECT v3_2 software as described by 
IBlanton et ail (|2003fJ) . 

We use an empirically-determined redshift- dependent 
colour separator of u — r = 2.1 + 4. 2z, where we use observer 
frame rather than rest-frame colours; within these luminos- 
ity bins, the fraction of red galaxies is 0.40 (L3), 0.52 (L4), 
0.64 (L5), and 0.80 (L6). 

3.2 SDSS spectroscopic LRGs 

Ano ther sample we use is t he DR4 spectroscopic LRG sam- 
ple (|Eisenstein et al]|200lf) . This sample has a fainter flux 
limit (r < 19) than the Main sample, and colour cuts to iso- 
late LRGs. We include these galaxies in the redshift range 
0.16 < z < 0.35, for which the sample is approximately 
volume-limited and includes 36 278 galaxies total. 

In order to study variation within this sample, we 
use cuts on several parameters. First, we construct lu- 
minosities using the r-band model magnitudes, and de- 
fine three luminosity subsamples as shown in Table [T] 
The absolute magnitude cuts are defined in terms of h = 
_ffo/(100 kms -1 Mpc -1 ) such that one can implement the 
cuts without specifying the value of Ho. The magnitudes 
have been corrected for extinction, and are k + e-correcte d 
to z = using the same templates as in IWake et al.1 |2006f) . 
We exclude galaxies lying inside the bright star mask. Ran- 
dom catalogs were generated taking into account the varia- 
tion of spectroscopic completeness with position. The ran- 
dom points were assigned absolute magnitudes and redshifts 
drawn from a distribution consistent with the real sample. 
Shape measurements were obtained via re-Gaussianization 
for 96 per cent of this sample (see i)3.4p : the remainder failed 
due to various problems, such as interpolated or saturated 
centers. 
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Table 1. Description of different galaxy samples and their subsamples used for the intrinsic alignment measurements. Included is a 
sample code and a brief description of the sample, the redshift range spanned, the absolute magnitude and colour cuts (including the 
redshift to which the absolute magnitude was k + e-corrected), the magnitude type (P for Petrosian, M for model), and the number of 
galaxies. 



Sample 


Sample description 


z range 


M r Cuts 




m r type 


N gal 


L3 


SDSS Main L3 (Section l3~ll 


0.02 < z < 0.12 


-20 M0- 1 <C - 


19 


P 


66 312 


LI 


SDSS Main L4 (Section 13.111 


0.02 < z < 0.19 


-21 < M r hl < - 


20 


P 


118 618 


L5 


SDSS Main L5 (Section E} 


0.02 < z < 0.29 


-22 < M' r hl < - 


21 


P 


73 041 


L6 


SDSS Main L6 (Section [37TJ 


0.03 < z < 0.35 


-23 < M\ )A < - 


22 


P 


7 937 


LRG 


SDSS LRG ( Section 


0.16 < z < 0.35 






M 


36 278 


LRG1 


cut by magnitude 




M°-° ^ -22.3 






16 068 


LRG2 


cut by magnitude 




-22.6 M°-° < - 


22.3 




13 019 


LRG3 


cut by magnitude 




M°-° < -22.6 






7191 


LRG.BG 


cut by environment 




BG only 






30 849 


LRG.non-BG 


cut by environment 




non-BG only 






5 429 


2SLAQ 


Full 2SLAQ LRG (Section l3~3)l 


0.4 < z < 0.8 






M 


7758 


2SLAQf 


cut by magnitude 




M° ° > -22 




M 


3 768 


2SLAQb 


cut by magnitude 




M°-° < -22 




M 


3 490 



Besides luminosities, we also use measures of local en- 
vironment to understand variation of intrinsic alignments 
within the sample. In particular, one question we may ask is 
whether intrinsic alignments are properties only of the cen- 
tral galaxy in a halo, or if satellite galaxies have intrinsic 
alignments as well. Unfortunately it is very difficult to actu- 
ally measure which galaxies are "central." A related quan- 
tity that we can measure is whether the galaxy is a brightest 
group/cluster galaxy (BG) or non-BG. To isolate BGs, we 
require that a given LRG be the brightest spectroscopic LRG 
within 2 Mpc projected separation and ±1200 km s _1 . 
This cut designates 83 per cent of the sample as BGs and 
17 per cent as non-BGs; the "BGs" are either in the field 
or host groups and clusters. One limitation to this cut is 
fiber collisions, for which the relevant radius is 55" (corre- 
sponding to ~ 200 h^ 1 kpc at the range of redshifts of the 
sample), which may cause us to mistakenly identify non-BGs 
as BGs. This is, however, only relevant on small scales. We 
emphasize that this distinction between BGs and non-BGs 
does not place any significant constraint on the "non-BG" 
galaxy luminosity, which may be quite close to the BG lu- 
minosity; the median luminosity gap is 0.3 magnitudes. Un- 
fortunately, the sample of non-BGs is too small to obtain 
meaningful results when splitting by luminosity gap. It is 
encouraging that our "non-BG" fraction of 0.17 agrees with 
the satellite fraction estim ated from halo mode ling of the 
LRG three-point function l|Kulkarni et al.l l2006). However, 
in interpreting the results, it must still be remembered that 
there could be cases where the BG is not actually the central 
galaxy (in some cases, such as merging clusters, there may 
not even be such a thing as a "central" galaxy), and that 
our selection is restricted to LRGs. In principle one could 
have increased the signal-to-noise ratio for the non-BGs by 
increasing the transverse separation cut to > 2/i _1 Mpc so 
as to obtain equal numbers in the BG and non-BG samples, 
but this would have had the undesired effect of cutting on 
large-scale environment (i.e. far beyond the virial radius), 
thereby including many central galaxies in the "non-BG" 
sample. 

Finally, in order to constrain redshift evolution more 



tightly, we have split each of the three LRG luminosity sam- 
ples at z = 0.27 to create a total of six spectroscopic LRG 
samples. 



3.3 2SLAQ LRGs 

The 2dF-SDSS LRG and QSO (2SLAQ) Survey has pro- 
duced samples of roughly 10 000 LRGs with z me d ~ 0.55 
with spectroscopic redshifts, and 10 000 faint quasars. For 
this paper, we use the 2SLAQ LRGs, which have similar 
cuts to the SDSS spectroscopic LRG sample, with a fainter 
apparent magnitude cut of rrt; (model) < 19.8 and with a 
slightly bluer rest-frame colour cut cj_ than in the SDSS 
LR G selection. 

ICannon et al. (2006) includes further details of the 
2SLAQ selection and observation; I Wake et al, I (|2006l) shows 
a first analysis of the LRG luminosity function evolution us- 
ing this sample in comparison with the SDSS spectroscopic 
LRG sample. While roughly 10 000 LRG redshifts have been 
measured, several of our criteria reduce the sample size: we 
require that they pass th e cuts to be in the LRG primary 
sample (|Wake et al.ll200?j . sample 8) rather than secondary 
sample (sample 9, which has far lower completeness); we 
eliminate repeat observations; we reject those in fields with 
significant incompleteness; we reject those with poor red- 
shift quality; we require that those with redshifts lie in the 
range 0.4 < z < 0.8; we require that they have shape mea- 
surements with sufficiently high resolution factor ( H3.4|l . The 
last of these cuts reduces the final sample from the canoni- 
cal 8 656 to 7 758 galaxies. Random catalogs with the same 
completeness as a function of angular and radial position as 
the real sample were used for the random points. 

Due to the faintness of this sample, we use model mag- 
nitudes. Thes e are k + e-corrected to z — using predictions 
derived from iBruzual fc Chariot! l|2003h ; the limitations of 
these models in describin g this sample are discussed more 
fully in I Wake et ail l|2006l ). 
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3.4 Ellipticity data 

In addition to a sample of galaxies, we also need their 
elli pticities. For this purpose , we use the measurements 
by iMandelbaum et al.l l|2005h . who obtained shapes for 
more than 30 million galaxies in the SDSS imaging data 
down to extinction-corrected magnitude r = 21.8, (i.e. far 
fainter than the spectroscopic limit of the SDSS). A mi- 
nor modification to the R eglens pipeline as described in 
IMandelbaum et all (|2006at ) was also used. 

The Reglens pipeline obtains galaxy images in the r 
and i filters from the SDSS "atlas images" (|Stoughton et all 
2002). The basic principle of shear measurement using these 
images is to fit a Gaussian profile with elliptical isophotes 
to the image, and define the components of the ellipticity 



(e+,e x ) = 



^^(cos2^,sin2^)), 



1 + (b/a) 



(8) 



where b/a is the axis ratio and (f> is the position angle of 
the major axis. The ellipticity is then an estimator for the 
shear, 



(7+,7x) 



2TZ 



<(e+,e x )), 



(9) 



where 1Z « 0.87 is called the "shear responsivity" and repre- 
sents the response of the ellipticity (Eq .[8l) to a small shear 
l|Kaiser et al.|[l995l ; lBernstein fc Jarvisl l2002). In practice, a 
number of corrections need to be applied to obtain the el- 
lipticity. The most important of these is the correction for 
the smea ring and circulariza t ion o f the galactic images by 
the PSF; IMandelbaum et al.1 (|2005l ) uses the PSF maps ob- 
taine d from stellar images by the PSP pipeline l|Lupton et al.l 
1200 j ). and c orrects for these using t he re-Gaussianization 
technique of iHirata fc Seliakl (|2003ah . which includes cor- 
rections for non-Gaussianity of both the galaxy profile and 
the PSF. A smaller correction is for the optical distortions 
in the telescope: ideally the mapping from the sky to the 
CCD is shape-preserving (conformal), but in reality this is 
not the case, resulting in a nonzero "camera shear." In the 
SDSS, this is a small effect (of order 0.1 per cent) which 
can be identified and removed using the astrometric solu- 
tion (|Pier et al.ll2003h . Finally, a variety of systematics tests 
are necessary to determine that the shear responsivity 1Z 
has in fact bee n determined corre c tly. W e refer the inter- 
ested reader to IMandelbaum et all (2005) for the details of 
these corrections and tests. 



4 METHODOLOGY 

The ba sic methodology in this pa per is very similar to 
that in IMandelbaum et al.l l|2006af) : we correlate a sam- 
ple of galaxies against some subset of itself, using the 
subset to trace intrinsic alignments and the full sample 
to trace the density field. We use the same set of esti- 
mat ors as in [Mandelbaum et al] (l2006aT ). generalizing the 
LS |Landv fc Szalavl 1 19931 ) estimator commonly used for 
galaxy-galaxy autocorrelations. 

The generalization of this estimator to the galaxy- 
intrinsic correlation is 



£ 9+ (r p ,II) = 



S+(D-R) _ S+D-S+R 



where S+D is the sum over all pairs with separations r p and 
II of the + component of shear: 



S+D : 



-CjIO 



211 



(11) 



RR 



RR 



(10) 



where e+(j\i) is the + component of the ellipticity of galaxy 
j measured relative to the direction to galaxy i, and TZ is the 
shear responsivity. S+R is defined by a similar equation. We 
emphasize that positive £ 9 + indicates a tendency to point 
towards overdensities of galaxies (i.e., radial alignment, the 
opposite of the convention in galaxy-galaxy lensing that pos- 
itive shear indicates tangential alignment). 

We also generalize this large-scale correlation function 
estimator to cross-correlations. In this case, we use one sam- 
ple to trace the intrinsic shear, and other to trace the density 
field. Thus in Eq. (|10[) . we find pairs of galaxies such that 
one is in the shear sample and the other in the density sam- 
ple, so that the S+ is determined from the former, and the 
D — R from the latter; the RR in the bottom is determined 
using one random point corresponding to the shear sample 
and the other corresponding to the density sample. 

This correlation function estimator is then integrated 
along the line of sight to form the projected intrinsic shear - 
density correlation function w g +(r p ). We model this function 
as a power-law, w g + — Ar p , where fits for A and a are done 
using the full jackknife covariance matrix. 

F or this paper, we use d the sa me two software pipelines 
as in Mandelb aum et al.l |200 6a) to compute the large- 
scale density-shear correlations. We independently general- 
ized both pipelines to allow different galaxy samples to be 
used to trace the intrinsic shear and the density field. In 
summary, the pipeline for which we present results com- 
putes the correlation functions over a 120 ft -1 Mpc (comov- 
ing) range along the line of sight (—60 < II < +60 h~ 1 Mpc) 
divided into 30 bins, then integrates over II. The range of 
transverse separations is from 0.3 to 60 hT x Mpc, in 10 log- 
arithmic bins. Covariance matrices were determined using 
a jackknife with 50 regions; as demonstrated there, results 
were relatively robust to the number of bins, with 100 re- 
gions giving approximately the same detection significance 
as 50 regions. For more d etails about these pipelines, see 
IMandelbaum et al.l (|2006al ). 



5 RESULTS 

In this section we describe the results of the galaxy density- 
shape correlation functions for the SDSS Main samples, the 
SDSS LRGs, and the 2SLAQ LRGs. The measured correla- 
tion functions and their uncertainties and correlation matri- 
ces are presented in Appendix |A1 



5.1 SDSS Main sample 

Here we present results of the measurements of the galaxy 
density-shape correlations w g +(r p ) using the Main sample 
split not only into luminosity bins, but also into colour sam- 
ples. There are 8 subsamples total, since we have 4 lumi- 
nosity bins and 2 colour bins. We attempt to address the 
issue of whether blue galaxies show any density-shape align- 
ment when correlated against the full L3-L6 (all colours), 
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w g+ for Main colour sub-samples 
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Figure 1. The GI correlation functions for the SDSS Main sub- 
samples, split into colour and luminosity bins. 
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Figure 2. Confidence contours for power-law fits to w g ^.(r p ) for 
Main sample galaxies. Contours are shown for various subsets of 
data labelled on the plots; in each case, la, 2a, and 3<r contours 
are shown. 

and to place constraints if there is no detection. The GI 
signal w g +(r p ) for each of the eight subsamples is shown 
in Figure [T] Figure [2] shows the confidence contours for 
fits to a power law, w + (r v ) = Ar p a as discussed in 
iMandelbaum et all l|2006al ). 

As shown, for both colour subsamples there is hardly 
any detection in L3 and L4, consistent with previous works. 
There is a hint of a signal in L3 for the red sample: when 
we fit the whole range of scales to an arbitrary power law 
there is no detection, however as we will see in Section [7] 



w g+ for spectroscopic LRG subsamples 




1 10 
r [Mpc/h] 



Figure 3. The density-shape correlation function w g +(r p ) from 
1-60 /i" 1 Mpc with the full spectroscopic LRG sample and various 
subsamples as labeled on the plot. Errors are la but are somewhat 
correlated on large scales. 

if we fit to large scales (r p > 4.7/i -1 Mpc) where the bias 
is expected to be roughly linear, and restrict to the power 
law a = —0.73 observed for the LRGs, there is a marginal 
(2.4<r) detection. For L5, the detection with the red sub- 
sample is robust whereas there is no detection with the blue 
subsample. For L6, the constraints are weak with the blue 
sample due to its small size, so while the magnitude of the 
alignments are consistent with the red subsample, they are 
also marginally consistent with zero. The rest-frame colour 
distribution of the L6 blue sample, and the distribution of 
Photo pipeline output frac_deV (a measure of the degree 
to which a galaxy profile is closer to an exponential or de 
Vaucouleur profile), suggest that this small L6 blue sample 
may contain galaxies that are on the edge of the blue vs. 
red galaxy distinction, which could explain this consistency 
of results. 

5.2 SDSS LRGs 

Here we present results of the measurements of density- 
shape correlations for the SDSS spectroscopic LRGs 
(the measurement is otherwise similar to that of 
IMandelbaum et al.ll2006ah . The plots of w g +(r p ) in Figure[3] 
are in the same form as in that paper, including la errors. 
Figure [4] shows the confidence contours for fits to a power- 
law, w g+ — Ar p . 

In the top panel of Figure[3] we show w g + (r p ) for the full 
spectroscopic LRG sample, and the BG and non-BG sub- 
samples. As shown, the full sample and the BG subsamples 
are robustly detected on all scales, out to 60 /i" 1 Mpc. The 
non-BG subsample has a significantly lower signal to noise 
due to its small size, but the amplitude appears roughly 
comparable to that of the BG subsamples. In Figure U the 
contours for the full LRG sample are shown in the upper 
left panel, and for BGs and non-BGs separately in the lower 
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Confidence contours for w g+ (r p ) = Ar£ 
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Figure 4. Confidence contours for power-law fits to w g +(r p ) for SDSS spectroscopic LRGs. Contours are shown for various subsets of 
data labelled on the plots; in each case, lc, 2cr, and 3cr contours are shown. 



left panel. As shown there, the constraints on the non-BG 
sample are indeed quite weak, but the detection is still ro- 
bust at slightly higher than the 3er level. The amplitude is 
higher and the power-law steeper than for the BG sample, 
but is still consistent with it at the la level. 

In the middle panel of Figure[3] we show w g +(r p ) for the 
spectroscopic LRG sample split at the median redshift. The 
purpose of this test is that, since this sample is roughly vol- 
ume limited, with the same luminosity-distribution in each 
of the halves, any evolution in the intrinsic alignment am- 
plitude should be due to redshift evolution. While the am- 
plitude appears higher for the higher redshift sample, the 
difference is not large; we must see the confidence contours 
on the fit parameters (which take into account the error 
correlations) before deciding if the difference is statistically 
significant, and will also later account for evolution of bias 
with redshift in Section [7] These contours appear in the up- 
per right panel of Figure [4] and indicate that the difference 
is, indeed, not statistically significant. 

Finally, the bottom panel of Figure [3] shows w g +(r p ) for 
the spectroscopic LRG sample split into luminosity subsam- 



ples. These results appear to have roughly the same power- 
law index, with a multiplicative difference in amplitude, and 
the confidence contours in the lower right panel of Figure [4] 
indicate that this is indeed the case. 



5.3 2SLAQ LRGs 

In this section, we present results of the measurement of 
the density-shape correlation using the 2SLAQ LRG sam- 
ple. Unfortunately, we cannot obtain meaningful results by 
comparing the full 2SLAQ sample results against the full 
spectroscopic LRG sample, because the samples do not cover 
the same range of absolute magnitudes. The 2SLAQ sample 
is, on average, fainter and bluer than the spectroscopic LRG 
sample. This difference is problematic since it will tend to 
give a lower intrinsic alignment amplitude for the 2SLAQ 
sample that we will need to separate from redshift evolution 
effects. 

To illustrate this point, Figure [S] shows a comparison of 
the k + e-corrected (to 2 = 0) model r-band absolute magni- 
tude distributions of the samples. First, the upper left panel, 
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Comparison: SDSS spectroscopic LRGs versus 2SLAQ LRGs 
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Figure 5. A comparison of the r-band absolute magnitudes and 
colours of the SDSS and 2SLAQ LRG samples. The upper left 
panel shows the absolute magnitudes as a function of redshift. 
The upper right shows the distribution of g — i rest frame colour, 
and the lower right shows the absolute magnitude histograms over 
the full samples. 



which shows scatter plots of magnitude M° versus redshift 
z. As shown, the spectroscopic LRG sample has roughly con- 
stant mean absolute magnitude as a function of redshift for 
the full redshift range used here. The 2SLAQ sample, on the 
other hand, shows a trend of being fainter at the low redshift 
end and brighter at the high redshift end. The mean absolute 
magnitude is roughly 0.2 magnitudes fainter for the 2SLAQ 
sample than for the SDSS spectroscopic LRG sample (this 
can also be seen in the bottom right panel, which shows his- 
tograms of the absolute magnitude values). This difference 
is a concern for our work, since in the previous subsection 
we showed that the intrinsic alignment amplitude has a very 
strong scaling with luminosity. 

One might wonder whether this difference that we ob- 
serve is due to s ome systematic, su ch as uncertainty in the 
k+e-corrections. IWake et alJ (|2006h explore the difference in 
more detail, and it does not seem possible to design a plau- 
sible model with different formation times, levels of star for- 
mation, or other differences that would reduce the difference 
between the samples to zero. They are inherently different. 

Because of this difference between the samples, rather 
than comparing intrinsic alignment amplitudes directly, we 
choose instead to fit to a model of the intrinsic alignment 
amplitude as a function of luminosity and redshift sepa- 
rately. This work can be aided by splitting the 2SLAQ sam- 
ple into bright and fainter subsets at M r — — 22 in addition 
to using the six redshift and luminosity samples of SDSS 
spectroscopic LRGs. Figure |§] shows r p w g +(r p ) for the two 
2SLAQ luminosity samples, and Figure [7] shows the confi- 
dence contours for a fits to a power-law for each subsample. 

Results for the luminosity and redshift-dependent fits 
will be presented in Section [7] 



Figure 6. The galaxy density-shape correlation function 
r p w g j r (r p ) from 1—60 h~ 1 Mpc with the luminosity subsamplcs 
used to trace the shapes and the full 2SLAQ sample used for the 
galaxy density g. Errors are lcr but are slightly correlated on large 
scales. 
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Figure 7. Confidence contours for power-law fits to u) 9 +(r p ). lcr, 
2er, and 3<r contours are shown for both luminosity subsamples. 

5.4 Systematics tests 

As in lMandelbaum et al.l l|2006al ). we have done several sys- 
tematics tests to ensure that the detections have not been 
contaminated by spurious instrumental or other effects. The 
first is the standard 45-degree test, whereby we rotate all 
source ellipticities by 45 degrees before computing the corre- 
lation functions (to obtain w gx instead of w g +). This rotated 
correlation function reverses sign under parity (i.e. the sign 
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of the correlation function is flipped depending on whether 
one rotates clockwise or counterclockwise) and therefore 
cannot appear unless there is a systematic, or galaxy for- 
mation violates parity invariance|j This test was done for 
the six subsamples of the SDSS spectroscopic LRGs (split- 
ting jointly into three luminosity and two redshift bins), for 
the "BG" and "non-BG" subsamples of the SDSS LRGs, 
and for the two luminosity subsamples of the 2SLAQ LRGs, 
giving ten w gx computations total. For each computed w gx 
we computed the x 2 f° r a fit to zero signal, and found the 
associated p- values (the probability of g etting a larger value 
of y 2 b y chance). We note that, as in iMandelbaum et al.l 
l|2006al ). the \ 2 values do not follow the usual \ 2 distribu- 
tion because of noise in the jackknife covariance matrices. 
The formalism to describe thi s effect and account for it with 
simulations was developed in iHirata et ail (2004). The ten 
p(> x 2 )-values computed taking into account the modified 
distribution of \ 2 ranged from 0.2-0.94; the lack of very 
low values indicates that we do not have significant B-mode 
contamination. 

Another test from IMandelbaum et al.l (l2006al ) is the in- 
tegration of the 3-dimensional correlations £ 9 + over large 
line-of-sight separations, 30 < |An| < 90 Mpc. A sig- 
nal that is of astrophysical origin should be dominated by 
the results from smaller |AII|, and hence should be consis- 
tent with zero for this range. This test will allow us to rule 
out, e.g., some optical effect causing an apparent alignment 
of galaxy images (which would not depend on the relative 
line-of-sight separation). For this test, in all cases, the sig- 
nal when integrating over large line-of-sight separations was 
consistent with zero, with p(> \ 2 ) values ranging from 0.07 
to 0.94. 

As infMandclb aum et al.l (|2006al ) , we return to the ques- 
tion of the size of jackknife subsamples. This was primar- 
ily a concern for the lower redshift subsamples, for which 
a jackknife region of a particular size corresponded to a 
small comoving transverse separation. Nonetheless, despite 
the higher average redshift of the samples in this paper, 
we revisit this issue briefly. A particular concern is for the 
2SLAQ sample which, while at higher redshift than any of 
the other samples, also covers a much smaller area overall, 
so independence of the jackknife regions may be a concern. 
For the samples used here, changing the number of jackknife 
samples by a factor of two did not change the size of the er- 
rors by more than 5-10 per cent, which suggests that we are 
using the jackknife in a regime where it has converged. 

A final concern we may have in comparing results from 
the SDSS spectroscopic LRG at z ~ 0.3 and the 2SLAQ 
LRGs at 2 ~ 0.55 is that the shape measurements, which 
are averaged over r and i bands, are for different sets of rest- 
frame wavelengths. In principle, these may then be domi- 
nated by different stellar populations within the same galaxy 
that have different intrinsic alignment properties. To test for 
this effect, we have calculated the analogy of the GI correla- 
tion function in Eq. (|10|l for SDSS LRGs using the difference 
between r- and i-band shears rather than their average as 



3 Note that for intrinsic-intrinsic ellpiticity correlations the ro- 
tated correlation function w x x is parity-allowed and so this would 
not be a good systematics test; rather one would use mj_|_ x , i.e. 
rotate the ellipticity of only one of the two galaxies in question. 



is usually done. Assuming that intrinsic alignments do not 
vary across this range of wavelengths, this quantity should 
then be consistent with zero. (We remind the reader that 
r band for the SDSS LRGs corresponds roughly to i band 
for the 2SLAQ LRGs when considering the rest-frame wave- 
lengths.) For the six SDSS LRG subsamples, these correla- 
tion functions were consistent with zero, with p(> x 2 ) values 
ranging from 0.29-0.81. 



6 BIAS OF DENSITY TRACERS 

To convert the observed w g + to the more relevant (for cosmic 
shear) ws+, we need to know the bias of the sample used 
as the density tracer. This needs to be done for the two 
LRG samples, and for the SDSS Main sample. Due to the 
much higher signal-to-noise in the LRG samples, we have 
measured the bias from the LRG projected autocorrelation 
function. The SDSS Main sample has a much lower signal-to- 
noise detection of GI alignments, and so we have constructed 
a crude bias estimate by combining previously published 
results. 



6.1 LRGs 

Besides computing w g +, our code computes w gg , the galaxy- 
galaxy autocorrelations. We have compared our results for 
the full SDSS s pectro scopi c LRG sample wit h those in 
lEisenstein et al.l i|2005l ) and IZehavi et ail (|2005T l and found 
agreement at the la level. Likewise, our resu lts for the 
2SLA Q sample are consistent with those in iRoss et"aL1 
(2006), though the central value of derived bias does not pre- 
cisely agree due to significantly different analysis method- 
ologies. We determine consistency by computing the best- 
fit powerlaw to w gg and comparing the \ 2 f° r °ur best- 
fit parameters ve rsus the \ 2 f° r the best-fit parameters in 
IRoss et alj (|2006l ). and find p(> Ax 2 ) = 0.2. 

The projected autocorrelations are computed for 0.3 < 
r p < 60 h^ 1 Mpc using 

W gg (r P ) = J U(r P ,n)<ffl, (12) 

where £ gg has been estimated using the LS estimator, 
(D - Rf DD ~ 2DR + RR 



(13) 



RR RR 
Integration along the II direction is carried out for —60 < 
II < +60 hT 1 Mpc. As for the w g + calculations, covariance 
matrices are determined using a jackknife with 50 regions. 

We did three types of fits to obtain the bias from w gg for 
each of the four samples (low-z SDSS LRGs, high- 2 SDSS 
LRGs, faint 2SLAQ LRGs, and bright 2SLAQ LRGs). The 
first method is to fit the bias b to the linear correlation 
function, 



w gg {r p ) = b 



2 J kdk 

~2^r 



Pun(k)W rp (k)+C. 



(14) 



The s econd method is to use the Q-model from ICole et al] 

mi), 



2 I kdk , .1 + Qk" 
w gg( r P ) = / -r— ilin(fc)- 



2?r 



1 + Ak 



W rp (k) + C, 



(15) 



where for a real-space power spectrum A = 1.7h 1 Mpc 
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|Cole et alj|2005l ): the \ 2 is minimized with respect to both 
fe and Q. The third method is to use Eq. (|14jl but with 
the nonlinear power spectrum instead of the linear. The 
window function for two-dimensional projection for an in- 
finitesimally thin range in radius is simply a Bessel func- 
tion: W Tp (k) — Jo(kr p ). For a finite range in r p , the window 
function is the weighted average (by area) over r p : 



W rp (k) 



!r P '2in 2nr p Jo(kr p ) dr p 
J^ p max 2nr p dr p 



2[Ji(fer p , n 



Jl (^^p,min)] 



nm) 



(16) 



The constant C accounts for the effect of the integral con- 
straint on the numerator of Eq. (|13[1 . There are in principle 
additional correc tions associated with the uncertainty in the 
denominator (see lHui fe Gazta naga 1999 for a thorough dis- 
cussion), which are relevant if the correlation function is of 
order unity on the scale of the survey region. This is not the 
case here so we have not included them. 

We used the "WMAP+all" ACDM cosmology 
(|Spergel et all 120071 ) for the bias determinations, which 
has n m = 0.262, h = 0.708, Q b = 0.0437, n s = 0.938, 
and erg = 0.751; in general for other values of erg, the bias 
scales as b oc a^ 1 . The projected correlation functions and 
best-fit Q-m odels are shown i n Figu re [8] The linear theory 
fit used the lEisenstein fc Hul (| 19981 ) transfer function; the 
two nonlinear fits differ only in their use o f t he nonlinear 
mapp ings of either iPeacock fc Doddsl l | 19961 ) or lSmith et al.l 
|2003l ). 

The biases of the LSS tracers are displayed in Tabled 
note that two values of the minimum r p were used, with 
the linear fit restricted to the largest scales. The b g values 
are consistent with each other and the \ 2 values are rea- 
sonable. The Q values obtained with the Q-model fits are 
61 ±48 (SDSS Iow-z), -24 ±49 (SDSS high-z), and 50 ±128 
(2SLAQ), i.e. Q is poorly constrained at these scales but 
is con sistent with other determi nations for similar galaxy 
types (|Padmanabhan et al.ll2006T ). 

As a systematics test, we also computed the integral 
/£(Vp,n)dn over tne range so < |n| < 120/i _1 Mpc and 
compute d b 2 a using this in place of w p (r p ). (For the fit we 
used the ISmith et al.l 120031 nonlinear mapping.) In the ab- 
sence of systematics, this should give a result close to zero 
since the galaxies considered are well separated along the 
line of sight, but could be nonzero if (for example) there are 
spurious angular fluctuations in the galaxy distribution. The 
procedure gives "fe^" values of ±0.07 ± 6.15, -0.06 ± 0.16, 
and ±0.17±0.34 (la) for the SDSS low-z, SDSS high-z, and 
2SLAQ samples, respectively. These are all consistent with 
zero and much less than the actual b 2 g values. 

For the rest of the paper we have used the b iases from 
the 7.5-47/1" 1 Mpc fit with the lSmith et alj l|2003l ) nonlinear 
power spectrum. 

6.2 SDSS Main sample 

The density tracer for the SDSS Main sample is the com- 
bined L3-6 sample. The bias of this tracer varies as a func- 
tion of redshift because the nearby part of the sample is 
dominated by galaxies with L ^ and the more distant 
part is dominated by galaxies with L > L*, which show 



Table 2. The bias of the various tracer samples for four fits. 
The error estimates arc la. The fittin g methods are "li n" (linear 
theory); "Q" (Q-m odel); "PD" jPeacock fc Doddslll996h ; and "S" 
(Smith et al.ll2003l) . 



Sample 


r p range 


Method 


bg 


X /dof 




Mpc 








SDSS low-2 


12-47 


lin 


1.84 ±0.15 


0.924/1 


(z) = 0.22 


7.5-47 


Q 


2.05 ± 0.12 


0.809/1 




7.5-47 


PD 


1.97 ±0.12 


1.476/2 




7.5-47 


S 


2.01 ±0.12 


1.428/2 


SDSS high-z 


12-47 


lin 


1.94 ±0.09 


0.001/1 


(z) = 0.31 


7.5-47 


Q 


1.97 ±0.08 


0.068/1 




7.5-47 


PD 


1.93 ±0.07 


1.405/2 




7.5-47 


S 


1.97 ±0.07 


1.257/2 


2SLAQ 


12-47 


lin 


1.98 ±0.25 


0.879/1 


(z) = 0.55 


7.5-47 


Q 


2.22 ±0.24 


0.798/1 




7.5-47 


PD 


2.10 ±0.20 


1.026/2 




7.5-47 


S 


2.13 ±0.20 


0.953/2 



stronger clustering. This should be taken into account when 
converting w g +(r p ) to wg + (r p ); in particular it results in 
an effective bias fe e ff that varies according to sample whose 
intrinsic alignments are being measured: the L3 GI corre- 
lation function measurement is dominated by nearby pairs 
of galaxies, hence 6 e ff is low, while for L6 fe e ff is high. A 
very crude way of estimating this effective bias is to take 
the pair-weighted average of the bias of the density tracer, 



/ bgri 



com „ com 



dV 



,com n com d y 



(17) 



where bs is the bias of the density tracer, and n^ om and 
n^ om are the comoving number densities of the ellipticity 
and density tracers, respectively. For the bias of the de nsity 
tracers, we use the results from iTegmark et~aH (|2004T ) and 
ISeliak et all i|2005l ), namely that ogfe = 0.764 (L3), 0.848 
(L4), 0.968 (L5), and 1.427 (L6). The resulting values of 
fe e / / for each of the ellipticity tracers are shown in Table [3] 
The biases of the L3-6 samples are generally determined to 
an uncertainty of several per cent; we do not explicitly prop- 
agate these uncertainties as they will contribute negligibly 
to the final error in ws+(r p ). 



7 POWER LAW FITS FOR LRGS 

In this section we present an empirical parameterized model 
for the large-scale GI correlation of LRGs, with power-law 
dependence on the galaxy's luminosity, redshift (technically 
1 ± z), and transverse separation. It should be noted that 
such a model is purely empirical, in the sense that it repro- 
duces the correct correlation between the GI amplitude and 
the specified galaxy properties. This is all that is needed 
to predict GI contamination for cosmic shear surveys; how- 
ever when considering theoretical interpretation of the re- 
sults presented here, it must be remembered that our fits 
represent correlations, not causal relationships. For exam- 
ple, we find that the GI amplitude scales as roughly L , 
but this does not mean that the alignment is "caused by" 
the galaxy's r-band luminosity - in practice both the GI 
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Figure 8. The projected g alaxy-galaxy correlation functions w gg (r p ) for the three samples and the best-fit nonlinear power spectrum 
(using the Smith ct al. 2003 mapping) to the 7.5— 47h~ 1 Mpc. Note that the error bars are correlated. 



Table 3. The effective bias fe e fj of the L3— 6 density tracer used in 
the GI correlation analysis. As explained in the text, b c g depends 
on the SDSS Main subsample used to trace the ellipticity field 
because the sample is not volume-limited. Values are normalized 
to WMAP (<rg = 0.751); for other values of erg, these values should 
be re-scaled as fe c g oc rjg 1 . The weighted effective redshift z e g at 
which the correlation is measured is defined by replacing 6,5 in 
Eq. JT7} with z. 



Ellipticity tracer 



ell 



Zcff 



L3.red 
L3.blue 
L4.red 
L4.blue 
L5.red 
L 5. blue 
L6.red 
L6.blue 



1.08 
1.08 
1.11 
1.12 
1.14 
1.16 
1.16 
1.19 



0.07 
0.07 
0.09 
0.09 
0.10 
0.12 
0.12 
0.13 



amplitude and the galaxy luminosity are both products of 
the process of galaxy formation and hence are correlated. 
In particular, different scalings might be obtained if we also 
fit dependence on e.g. colour, velocity dispersion, surface 
brightness, etc., since these are correlated with luminosity 
and redshift. 

The simplest model including luminosity, redshift, and 
scale dependence is 





ws+(r p ) = A 



1 + z 



(18) 



10' 



0.4(-22-M r ) 



.T pivot / / ivot 

Here there are 4 parameters {Ao,a, j3, 7}, and L is the 
galaxy luminosity. We use the r-band luminosity k + e- 
corrected to 2 = 0.0. The normalization Lq corresponds to 
absolute magnitude —22, i.e. 

L_ 

T 

For reference, typical k + e-corrections for a completely 
passive z = 0.3 LRG are (1.51,1.00,0.15,-0.05,-0.11) 
in the five bands, and for a z — 0.55 LRG these are 
(2.73,1.59,0.77,0.05,-0.15). If the LRG has some con- 
stant low level of star formation, the k + e-corrections 
are (0.86,0.92,0.15,-0.04,-0.10) from z = 0.3 and 



(19) 



(1.04, 1.40, 0.73, 0.05, -0.13) from z = 0.55. In practice, k + e 
corrections are obtained by considering the observed value 
of g — i for each galaxy, and doing linear interpolation us- 
ing its relationship to the predicted observer-frame g — i 
values for the completely passive and the passive plus star 
forming models at that redshift (at z — 0.3, the models pre- 
dict g — i = 2.30 and 2.14 for the two models, respectively; 
at z = 0.55, they predict g — i — 2.79 and 2.53). The pivot 
points selected here are z p i V ot = 0.3 and r p i vot = 20ft. -1 Mpc. 



We have performed a least-squares fit of this ws+(r p ) 
to eight LRG samples: the SDSS LRGs split into six sam- 
ples (the 3 bins in luminosity and 2 in redshift), and the 
2SLAQ LRGs split into a faint and bright sample (sepa- 
rated at M r — —22.0). We have converted the observable 
w 9+( r p) to wg+(r p ) by dividing by the bias of the density 
tracer. For details on the bias determinations, see Section HI] 
Note that the bias determination assumed as = 0.751; if one 
wishes to have a constraint on u>s+ for other values of as it 
is necessary to multiply the measurement of Aq in Table [3] 
by (7 8 /0.751. 



The relation ws+(r p ) — w g +(r p ) /b g is the "obvious" 
way to do the conversion from correlation functions involv- 
ing galaxies to those involving matter, but it is worth con- 
sidering the specific assumptions this makes. It is valid pro- 
vided that (i) galaxies are locally linearly biased, i.e. the 
probability distribution for g(r) conditioned on a particu- 
lar realization of the density field depends only on S(r) and 
not on S at other locations, and that the dependence of the 
mean galaxy density is linear, (g{r))\s = b g 5(r); and (ii) the 
intrinsic shear of a galaxy depends only on the surrounding 
density field (or equivalently on the tidal field, which con- 
tains the same information) and not on the placement of 
the galaxies, P(~f I \8) — P(7 7 |5, g) (i.e. specifying the po- 
sitions of other galaxies provides no additional information 
not in the density field). Note that the dependence of in- 
trinsic shear on the density field is allowed to be nonlocal. 
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Given these two assumptions we may write, for r^O, 
£ 9+ (r) = (7^(0) ff (r)) = ([l + 5(0)hi(0) fl (r)) 

r V5 Vg P(5, fl )(7^(0))| 5 , g [l +g(0)]g{T) 
V5 P(*)(7i(0)>| 4 <[l + ff(0)]fl(r))| 4 



2>«5 P(5)(7^(0))| 5 (l+ fl (0)) a (ff(r)>|, 
2>«5 P(S)( 7 ;(0)[l + j(0)]) f M(r) 



b 9 <7+(0)[l + .g(0)]5(r)) = & 9 6+W. 



(20) 



Here the first two lines only involve defintions; the third line 
uses assumption (ii) to remove the conditioning of {7+ (0)) 
on g; the fourth line makes use of assumption (i) that the 
galaxy density depends locally on the matter density to split 
the expectation value ([1 + g(0)]g(r))\s; the fifth line uses 
the linearity of the biasing from assumption (i) to introduce 
b g , and uses assumption (ii) to combine two of the expecta- 
tion values into (7 + (0)[l + g(0)]}s; and the last line again 
only use definitions. The projected (2-dimensional) relation 
w s+( r p) = w g+{ r p)/b g then follows. (Here J T>5 represents 
a functional integral over realizations of the density field.) 
Gaussianity is not assumed, however Gaussianity combined 
with scale-independent bias (defined by b g = \Z£gg/£ss) 
and unit stochasticity r g = 1 would imply the validity of 
assumption (i) since they completely specify the two-point 
behavior of the galaxy field. We also note that (i) implies 
b g ^constant and r g = 1. Of our assumptions, (i) is consis- 
tent with (though not uniquely proven by) simulation results 
suggesting b q =constan t and r g = 1 down to ~ 5ft -1 Mpc 
l|Tasitsiomi et al]|2004 ); at smaller scales there is a hint of 
decrease in the bias. The status of assumption (ii) is un- 
known, however since galaxy alignments are expected to be 
determined by tidal fields it is physically reasonable. In the 
future it would be desirable to do a more detailed analy- 
sis, perhaps fitting the w g+ (r p ) results directly to simula- 
tions; for this reason the correlation functions are given in 
Appendix [S] We have also not attempted to use the con- 
version ws+(r p ) — w g +(r p )/bg at smaller scales than the 
4.7-7.5/i _ Mpc bin, although these inner bins are provided 
in the data tables. 

We perform several fits to Eq. (|18p . In each case, since 
the galaxies used to trace the ellipticity field do not all have 
the same luminosity or redshift, we actually fit 



w s+ (r p ) = A - 



(l + ^pivot) 7 ' 



(21) 



pivot -^0 

where the luminosity and redshift averages are taken over 
the galaxies in the sample, and the r p average is weighted 
by area (analogously to Eq. I16|l . In principle one should also 
weight the redshift average by the comoving number density 
of the density tracer (full SDSS LRG or 2SLAQ sample) but 
in practice the comoving number density is nearly constant 
for the former and the statistical errors are very large for the 
latter, so we did not implement a correction for this effect. 

The fits were done including the correlations between 
different radial bins for the same subsample, but not in- 
cluding the correlations between different subsamples. The 
correlations between different radial bins of w g +(r p ) for the 
same subsample are clearly seen in any of the covariance 



matrices and are expected because galaxies are clustered. 
They must be included to get meaningful results. However, 
we do not expect significant correlations between w g +(r p ) in 
different subsamples, so long as the intrinsic alignments are 
sufficiently weak. In order to test this, we computed the cor- 
relation coefficients piA,jB by the jackknife procedure, where 
we have used the indices to denote radial bins and 

A,B, ... to denote LRG subsamples. We consider the 6 ra- 
dial bins used for fits here (4.7 < r p < 60/t~ 1 Mpc) and the 6 
SDSS LRG subsamples. This provides 315 cross-correlation 
coefficients with A 7^ B. In each case we then calculate the 
Fisher z-coefficient, defined by piA.jB ~ tanh ZiA.jB- These 
315 coefficients have a sample mean of +0.019 ± 0.008 and 
a standard deviation of 0.14. For comparison, with 50 jack- 
knife regions of the same size and Gaussian errors, we should 
have a standard deviation of ~ 1/^50 — 3 = 0.15, and the 
mean value of the {-ZiAjs} should be zero if the samples 
truly are independent. In fact we observe the correct stan- 
dard deviation and a hint of correlation between bins at only 
the ~ 2 per cent level. For the two 2SLAQ LRG bins there 
are 21 ^-coefficients with a sample mean of +0.016 ± 0.046 
and a standard deviation of 0.21. Therefore we have not 
included correlations between different samples. 

The fit results are shown in Table [4] We have included 
both SDSS-only and SDSS+2SLAQ constraints; it is clear 
from the fit that the SDSS LRG sample provides essentially 
all of the constraint on a and /3, with the addition of 2SLAQ 
providing significant new information about the redshift evo- 
lution parameter, 7. This is not surprising since the 2SLAQ 
sample is small but adds a large baseline in redshift. Note 
that the errorbars on Ao are rather non-Gaussian, so the 
apparent ~ 6a effect is actually stronger than that (see the 
Ax 2 values). 

The error bars in the table require some explanation. 
The usual way of constructing error bars on a single pa- 
rameter is to change its value until the \ 2 (minimized with 
respect to the other parameters) increases by an amount 
Ax 2 = 3.84 (the 95th percentile of the \ 2 distribution with 
1 degree of freedom). The problem with this procedure is 
that if the error bars are determined by resampling (jack- 
knife or bootstrap), the resulting covariance matrix C of 
the w g +(r p ) values is noisy (but still unbiased). This leads 
to a systematic overestim a te of C -1 l|Hirata fc Seliakl |2004| ; 
IMandelbaum et al.ll2006al ; lHartlap et al.ll2007h . In particu- 
lar, if one has a set of parameters {p 1 } being fit, the shift A\ 2 
(i.e. the x 2 using the true value of parameter i minimized 
over the other parameters, minus the x 2 minimized over all 
parameters, and using the estimated covariance matrix) is 
not distributed as a \ 2 with 1 degree of freedom, even for 
Gaussian-di stributed errors. A solu tion to this problem was 
propos ed bv lHirata fc Seliakl (|2004T ) and IMandelbaum et al.l 
l|2006al ) who showed that if the data have Gaussian errors, 
the covariance matrix is obtained from Gaussian data, and 
the model is linear in the parameters, then the distribution 
of Ax 2 does not depend on the covariance matrix C and 
hence can be obtained from a Monte Carlo simulation. Un- 
fortunately their procedure does not apply directly to our 
problem because we only estimated the covariance matrix 
entries between w g +(r p ) values for the same subsample - the 
cross-entries between different subsamples are zero. We have 
circumvented this problem by doing a Monte Carlo simula- 
tion as follows. First, we take the jackknife covariance matrix 
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C to be the fiducial model and simulate th e re-sampling by 
the pr ocedure described in Appendix D of iHirata fc Seliakl 
(2004) with M = 50 regions to get a re-sampled matrix 
C. We then set the cross-entries between different subsam- 
ples to zero. We also generate a simulated set of estimated 
{w g +(r p )} according to a Gaussian distribution with mean 
given by the best-fit power law model and covariance C. 
From this simulated data and simulated covariance matrix, 
one can then compute a x 2 surface and find Axi for each of 
the four parameters {p l }f =0 . Repeating this procedure 1024 
times allows us to compute the A\i distribution for each 
parameter. The 95th percentile of this distribution is shown 
in Table [4] for each fitting region and each parameter, and 
is used to compute the 95 per cent confidence limits on the 
parameters. 

The question naturally arises whether the fits to the 
LRG signal in Table [4] apply to the GI correlation for the 
red galaxies fainter than the LRG1 and 2SLAQ-faint sam- 
ples. Certainly in making a model for the GI correlations 
for cosmic shear purposes it is necessary to have a model for 
these fainter objects, which are after all much more numer- 
ous. The best way to assess this is to take our LRG models 
and compare them with the measurements for the L3, L4, 
L5, and L6 samples of red Main galaxies. It is seen from 
Figure [2] that our best-fit value of a = —0.73 is consistent 
with the slope in all four cases, so it remains to test the am- 
plitude. We do this by taking the GI measurements for each 
of these samples at r p > 4.7/i _1 Mpc and fitting to them a 
power- law 



Ws+(r p ) = raj + (r pivo t) 



(22) 



where a is fixed to —0.73 and the normalization ^5+ (^"pivot) 
is allowed to vary to minimize the x 2 ■ The resulting ampli- 
tudes are shown in Table [S] The L5 and L6 samples actually 
include some of the SDSS LRGs, so only the comparisons 
for L3 and L4 represent an independent test of the power- 
law model. There are detections at > 2a for the L3, L5, and 
L6 red samples, although the significance is far greater for 
L6 than for the others. We have also displayed, in the last 
column, the amplitude predicted by the best LRG fit, after 
converting the Petrosian magnitudes measured for the SDSS 
Main galaxies to model magnitudes by dividing by 0.8 (the 
fraction of the flux captu red by the Petrosia n method for 
a typical elliptical galaxy; iBlanton et alj|200ll ). This ampli- 
tude is consistent with all of the samples at the 95 per cent 
confidence level. 



8 COMPARISON TO THEORY AND 
SIMULATIONS 

In this section we compare our results to analytical and 
simulation-based models of intrinsic alignments. 



8.1 Analytical models 

It is worth comparin g the results of t his an alysis to the- 
oretical expectations. IHirata fc Seliakl <|2004h presented an 
analytic "linear passive" model for red galaxies, which as- 
sumes that a galaxy's ellipticity is a linear function of the 



Table 5. The GI correlation amplitude measured from the L3-L6 
red Main samples of galaxies, and that predicted by the best LRG 
fit (to SDSS+2SLAQ, r p > 4.7ri _1 Mpc data). The pivot radius 
is r pivot = 20h~ 1 Mpc. Note that the L5 and L6 samples have 
some overlaps with the LRG sample, so that for these samples 
the two columns are not independent. Errors shown are 95 per 
cent confidence limits, which for 50 jackknife regions, 6 radial 
bins, and 1 parameter (the amplitude) being fit corresponds to 
Ax 2 = 5.02. 



Sample 


Measured (95%CL) 


) (h- 1 Mpc) 
Predicted from LRGs 


L3.red 


+0.035 ±0.032 


+0.004 


L4.red 


+0.013 ±0.019 


+0.013 


L5.red 


+0.024 ± 0.022 


+0.044 


L6.red 


+0.144 + 0.078 


+0.132 



local tidal quadrupole when the galaxy forms and then re- 
mains unchanged. In this case we expect that on linear 
scales we should have P s *i (k, z) oc k n >> D(z), where nf is 
the scalar spectral index corrected for the transfer function, 
i.e. nf = n s + 2d [In T(k)]/d In k, and D(z) is the growth 
factor. The redshift dependence is simply the growth fac- 
tor because S increases but 7 does not, and the scale de- 
pendence is simply that of the underlying power spectrum. 
[This should be equivalent to Eq. (18) of IHirata fc Seliakl 
(2004), since there p oc (1 + zf, D(z) = (1 + z)D(z), and 
Ps in (k,z) oc k< !! D 2 {z). There is s till a factor of (1 + z) 2 
that is missing from Eq. (14-18) of IHirata fc Seliakl (|2004f ) 
because of the conversion from comoving to physical scale in 
the potential-density relation.] Transforming to real space, 
this predicts: 

n eff , dlnD(z) 
a = -2-n s and 7= d]n(1 + z) - (23) 

[There is a 2 instead of a 3 in a because ws+(r p ) is a pro- 
jected quantity and hence is obtained from P s ^7 (k) by a 
two-dimensional Fourier transform.] For a ACDM cosmol- 
ogy, this predicts that 7 should rise from — 1 at high redshift 
(matter domination) to in the far future (A domination). 
Across the range of redshifts considered here, 0.2-0.7, we 
expect 7 ~ —0.7. The prediction for a depends somewhat 
on scale since nf is not constant. We can find the value of a 
relevant to our observations by taking Ps^i{k) oc Pj lnear (fc) 
and using the Hankel transform 

w s+ (r p ) = -i- j P s ^(k)J 2 (kr p )kdk (24) 

to get u>s+(r p ). Measured across the range from (11.9- 
60)/i _1 Mpc (i.e. using the largest-scale cutoff from Ta- 
ble [4}, we find a — —0.65 for the fiducial cosmology. In- 
cluding smaller scales leads to an increase in a because the 
power spectrum curves downward, d 2 [lnP(fc)]/d(lnfc) 2 < 0: 
we have a = —0.52 over the range (7.5-60)/i _1 Mpc, and 
a = —0.41 over the range (4.7-60)fe _1 Mpc. The linear 
model does not give a prediction for /3, which requires an 
understanding of how the proportionality constant between 
ellipticity and tidal field relates to the luminosity. 

Our results for red galaxies are entirely consistent with 
the prediction 7 = 7 pa ssivc ~ —0.7, although given the large 
error bars in Table [4] this is not a particularly impressive 
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Table 4. The best-fit parameters to Eq. I|18| l using SDSS and 2SLAQ LRGs. Error bars are 95 per cent confidence limits. The first set of 
fit parameters uses the most conservative cut on transverse separation. The Ax 2 value in the second-to-last column is the improvement 
in x 2 relative to no GI correlation, and the Adof value indicates the number of parameters in the fit. The amplitude Ao is given in units 
of O.Ol/i -1 Mpc. The A% 2 values in the last column indicate the degradation of \ 2 used to compute the 95 per cent confidence limits on 
Aq, a, P, and 7, respectively. These are greater than 3.84 for 1 parameter because of noise in the covariance matrix (see text). 





Fit reg 


Ion 


Ao/CO.Olh- 1 Mpc) 


a 


p 


7 


X 2 /dof 


A x 2 /Adof 


AX? 












Fits to SDSS+ 


2SLAQ 








r p 
r p 
r p 


> 11.9/1 

> 7.5/i- 

> 4.7/i- 


_1 Mpc 
-1 Mpc 
-1 Mpc 


+6-0i 2 2 | 

+6.4± 2 .;i 
+5.9I 2 ;3 


-0.88±g;|i 


,1 ci +0.73 

-t-i.oi_ 69 

+1 41 +0.66 
+ l-48_ g3 


-1 oo+ 2 - 40 

-0.27jl|;|| 
-0.56I 2 ;? 2 


33.3/28 
42.8/36 
54.9/44 


171.8/4 
215.8/4 
219.2/4 


4.38, 4.74, 4.16, 4.47 
5.10, 4.89, 5.07, 5.02 
5.04, 5.08, 4.92, 5.57 


Fits to SDSS only 


r p 
r p 
r p 


> 11.9/1 

> 7.5/i- 

> 4.7/1" 


"iMpc 
- 1 Mpc 
-1 Mpc 


+7 4+ 2 ' 9 
+ '•^-2.4 

+6.6l 2 -j 


-o.95t° ;i 

-0.88t°; 24 

-o.74i8;ii 


+1.431^7? 

+i.3i±g : |S 

+1.44lg;g 


+1.941H 
4_2 qq+4-52 

4-1 81 + 4 - 52 


21.3/20 
27.9/26 
34.0/32 


173.9/4 
208.7/4 
213.3/4 


4.49, 4.67, 4.27, 4.27 
4.66, 5.13, 4.91, 4.66 
4.88, 5.10, 4.66, 5.00 



accomplishment. From a practical perspective, we have at 
least set an upper limit on 7, which enables upper bounds 
on contamination to be placed in cosmic shear investiga- 
tions. From a theoretical perspective, we are still unable to 
answer perhaps the most interesting question: is 7 larger or 
smaller than the passive evolution prediction? A value of 
7 > 7passiv C would imply that the intrinsic alignments of red 
galaxies were greater in the past, and were being reduced, 
perhaps due to relaxation, merger, or figure rotation pro- 
cesses that destroy pre-existing correlations. Conversely, a 
value of 7 < 7passiv C would imply that some process was 
causing the LRGs to align themselves with the large-scale 
density field, even at low redshift where LRGs are generally 
believed to be passively evolving. 

By contrast, the linear theory predictions for a are only 
in marginal agreement with observations, with the data giv- 
ing a smaller (more negative) a. The discrepancy is not sig- 
nificant if one uses the range (11.9-60)/i _1 Mpc; however one 
finds a discrepancy of Ax 2 = 9.02 (p — 0.005 using the dis- 
tribution from the Monte Carlo simulations) if one includes 
data down to 7.5/i _1 Mpc or Ax 2 = 13.11 (p < 10~ 3 ) down 
to 4.7/i -1 Mpc. The direction of this disagreement is what 
one would expect if nonlinear clustering on small scales en- 
hanced wg + (r p ), since this would tilt a to more negative 
values than predicted by linear theory. Indeed i f one substi- 
tutes the nonlinear power spectrum P s al (k) of ISmith et al.l 
( 2003) into Eq. (f2"4")l then the predicted values of a are —0.72, 
-0.68, and -0.67 for r p , min = 11.9, 7.5, and 4.7/i _1 Mpc re- 
spectively. These are in very good agreement la) with 
the measured slopes, however the theoretical justification for 
believing Pg~i (k) to trace the matter power spectrum in the 
nonlinear regime is dubious. 

The data on GI correlation for blue galaxies are too 
noisy to constrain their shape, so we have not attempted 
to compare this to theoretical models. We also note that 
the main theory for ali gnments of these galaxies, namely 
the tidal to rque theory (|Peeblesl 1 19691 ; iDoroshkevichl Il970l ; 
IWhitelll984r i. generically predicts zero GI correlation up to 
second order in pert urbation theory (|Hui fc Zhand |2002| ; 
iHirata fc Sehakll2004h . Thus more detailed theoretical cal- 
culations will be necessary to predict the scale dependence 
of ws+(r p ) for blue galaxies. 



8.2 Simulations 

The GI contamination to t he lensing signa l has b een in- 
vestigated in simulations by iHevmans et alj (|2006bj) . They 
used ACDM Af-body simul ations populated using the 
ICoorav fc~ Milosavlicvic (2005) conditional luminosity func- 
tion and several different models for determining the galaxy 
(as opposed to halo) ellipticity. They then computed the 
GI correlations using various models: an "elliptical" model 
for which galaxies were assigned the same ellipticity magni- 
tude and direction as the parent halos; a "spiral" model for 
which they were considered to be disks with some random 
misalignment with the parent halo angular momentum vec- 
tor; and a "mix" model with a mixture of the above. They 
concluded that only the mix model was consist ent with the 
measurements from Mandelb aum et al.l (|2006al ). 

Here it is possible to do a more detailed comparison. 
In Appendix [B] we derive the conversion from the intrin- 
sic shear-lensin g shear correlation functions measured by 
IHevmans et~aH (|2006b) to ws + (r p ). We have applied this 
conversion to the elliptical and mix models, and plotted 
these in Figure [9] Overplotted are the ws+{r p ) data for the 
LRGs (median magnitude ~ 2 magnitudes brighter than L*) 
and the L4 blue galaxies (i.e. ~ L+). Th e LRG data are con- 
sistent with the IHevmans et al.l l|2006bl ) "elliptical" model, 
which is the one with the strongest GI correlation. This may 
seem surprising since the "elliptical" model was designed to 
be a maximal estimate of GI, in the sense that the galax- 
ies were assumed to trace the ellipticity of their host haloes 
perfectly - one would expect that in the real Universe there 
would be some misalignment. However the "elliptical" model 
is also averaged over a range of halo masses. Since LRGs 
typically occupy the most massive haloes, and simulations 
have suggested that the d ensity-halo ellipticity correlation 
is strong er at high masses dHopkins et al.ll2005h . it is possi- 
ble that IHevmans et alj (|2006bl ) would have found a much 
stronger GI correlation if they had only considered LRGs. 
Simulation results broken into luminosity bins should thus 
be a priority: until they are availabl e we cannot say whether 
our LRG measurements confirm the IHevmans et all l|2006bl ) 
"elliptical" model, or if there is substantial misalignment in 
the real Universe that f ortuitously results in agreement with 
IHevmans et~ai1 |2006b) because of their lower typical halo 



© 0000 RAS, MNRAS 000, 000-000 



16 Hirata et al. 



GI correlations 





0.35 - 




0.3 - 




0.25 


Mpc) 


0.2 - 


x: 


0.15 - 






± 


0.1 - 


s 






0.05 




- 




-0.05 - 



L4 blue 
LRG 



H06 elliptical 



H06 mix 



3 5 8 10 20 30 50 

r p (If 1 Mpc) 

Figure 9. The GI correlation functions ui<5+(r p ) for the full SDSS 
LRG sample and the L4 blue sample, with the galaxy-to-density 
conversion using biases of 1.99 and 1.12 respectively. The lines are 
the curves from lHevmans et al. (2006b) at z = 0.27, interpreted 
using the procedure in Appendix [B] The median rcdshift of the 
LRGs 

is z med — 0.27; the L4 blue galaxies have z c ff — 0.12, 
although we expect that evolution between 0.12 and 0.27 would 
make a small effect. 



mass. In any case we note that LRGs show the strongest 
alignment signal of all galaxy types considered and make up 
only a small fraction of any flux-limited sample of galax- 
ies, so for a realistic weak lensing survey the GI contamina- 
tion would be lower. The more abundant L* blue galaxies 
are consistent ( within 2a) of either zero alignment or the 
iHevmans et"al] l|2006bD ' 'mix" model, and it is evident that 
more data will be needed to distinguish these possibilities. 

We are not aware of any simulation results that explore 
the luminosity dependence of the GI corre lation. Studies 
of clu ster alignments in iV-body simulations ((Hopkins et al.l 
l2005h find a preferential alignment of the major axis of the 
cluster with the direction to neighbouring clusters with a 
strength that increases with halo mass. This is qualitatively 
consistent with what we observe for the alignments of LRGs, 
which increase with luminosity, but a quantitative compari- 
son would require one to populate the haloes with galaxies, 
assign ellipticities, and convert alignment angles to density- 
shape correlation functions; each of these (but especially the 
latter two) is a significant source of uncertainty. 



9 ESTIMATES OF CONTAMINATION 

In principle, with knowledge of the GI correlation scaling 
with luminosity, colour, redshift, and transverse separation, 
combined with knowledge of the joint luminosity-colour- 
redshift distribution for a particular survey, we can predict 
the contribution of GI contamination to the measured cos- 
mic shear power spectrum. In principle, this contamination 
then carries over to an underestimate of as or, if the evolu- 
tion of the amplitude is measured to constrain the equation 
of state of dark energy, there will be errors in the measure- 
ment of wo and w a . In future work, we will attempt to quan- 
tify more precisely the effects on the dark energy parameter 
estimates. For now, we merely present the fractional contam- 
ination of the cosmic shear power spectrum measurement for 



model surveys with realistic redshift and luminosity distri- 
butions, and present a prescription for marginalizing over 
GI uncertainties in measurements of a%. 



9.1 Contamination models 

To make a prediction of contamination to cosmic shear re- 
sults, we start with a flux limit, and must assume our model 
fits for intrinsic alignment amplitude as a function of sepa- 
ration, redshift, and luminosity; a luminosity function as a 
function of spectral type; the underlying cosmology (to get 
the distance modulus); and k + e-corrections. From this, we 
can predict the redshift distribution and the magnitude of 
the GI correlation contribution to the measured cosmic shear 
power spectrum. In principle, there is also a correction asso- 
ciated with the types of magnitudes used, e.g. Petrosian vs. 
model. For example, the Petrosian system typically misses 
~ 20 per cent of the fl ux for well-resolved elliptical galax- 
ies l|Blanton et al.|[200lh . Since the GI signal for red galaxies 
scales as roughly L 1 ' 5 , this would translate into roughly a 30 
per cent error in the GI amplitude, which is small compared 
to the redshift extrapolation uncertainty in our models. To 
the extent that it matters, users of these models should note 
that our contamination estimate for red galaxies is fit to the 
LRGs, for which we used model magnitudes; for blue galax- 
ies the constraints on GI are so weak that changes of this 
order are unimportant. 

For the r-band luminosity function as a function of 
spectral type (including luminosity evolution), and for k- 
corrections, w e rely on the re sults from COMBO-17 for 
0.2 < z < 1.2 i|Wolf et al.ll2003h . For a typical cosmic shear 
survey that is dominated by galaxies with 0.8 < z < 1.2, 
the galaxies that are most important for the intrinsic align- 
ments are those at the lower end of the redshift range, where 
our assumptions about luminosity evolution, etc. are most 
likely to be valid. For the purposes of this work, the spec- 
tral types 1 and 2 in that paper are considered "red," types 
3 and 4 are considered "blue." The reason for this distinc- 
tion is that at z = 0, integration of the templates to obtain 
observed AB u — r colours (which were used in this paper 
for colour separation) suggests that all type 1 and nearly all 
type 2 galaxies would have been classified as red, whereas 
types 3 and 4 would have been classified as blue. At redshift 
0.1-0.25, the classification of the type 2 template as red ver- 
sus blue is no longer clear, but this uncertainty is in part due 
to certain features of the spectrum below 300 nm redshifting 
into the -u-band. These features have changed significantly 
in updated versions of the sp ectra used b y the COMBO-17 
team, since the templates of IWolf et al.l (|2003l ) were found 
to be inadequate for full galaxy classificatkirQ so we define 
our correspondance between SDSS and COMBO-17 types 
using the z = colour (which is not susceptible to influ- 
ence by these uncertainties in the A < 300 nm portion of the 
spectrum). With this classification scheme, roughly 25 per 
cent of the galaxies in a flux-limited survey to R = 24 are 
classified as red. 

The intrinsic alignment signal for the red galaxies is 
very well constrained at z ~ 0.3, with the major uncertainty 
being the scaling of the signal with redshift. The situation 
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is very different for the blue galaxies for which we have no 
detection, and for which we only have a low-redshift con- 
straint: the GI correlation for these objects may be near 
present upper limits, or alternatively could be negligible. In 
light of these uncertainties, we have defined four different 
models, all based on the fits at r v > 4.7/i _1 Mpc: 

1. In the "pessimistic" model (A), we assume that the red- 
shift scaling of GI for the red galaxies is the 95 per cent con- 
fidence upper limit from the SDSS+2SLAQ fit with r p > 4.7 

Mpc, 7 = +1.47. The GI values for the blue galaxies are 
taken to have the same radial scaling a and redshift scal- 
ing 7 as for the red galaxies, with the amplitude taken to 
be the 95 per cent confidence upper limit from SDSS Main. 
The use of the same radial scaling a is motivated by the 
expectation that the alignment on large scales should trace 
the tidal field, which has the same spatial dependence for 
all types of galaxies. A single amplitude is fit to the L3-5 
blue galaxies to improve statistics since there is no sign of a 
detection in any of these bins. A separate amplitude is used 
for the L6 blue galaxies as these have a marginal detection 
when a is fixed. 

2. In the "central" model (B), we use the best- fit values 
of amplitudes and power law indices a, f3, 7 from Table [3] 
using both SDSS and 2SLAQ and including radii down to 
A.7h~ l Mpc. For the blue galaxies we use the same value 
of a, as explained above, but fix 7 to the passive value, 
7 P assive = —0.7 and use the best-fit rather than worst-case 
value for the amplitude. The L3-5 galaxies are combined, 
just as for model B. 

3. In the "optimistic" model (C), we treat the red galaxies 
and the L6 blue galaxies just as for model B. We assume the 
L3-5 blue galaxies have no GI signal, since none is required 
by the data. (We do not view the detection in the L6 blue 
galaxies as robust since the significance is < 2.5cr and only 
appears in one of the four blue bins.) 

4. In the "very optimistic" model (D), we take the mini- 
mum value of 7 allowed by SDSS+2SLAQ at the 2a level 
(7 = —2.81), and use the values of the other parameters 
that minimize the \ 2 constrained to fixed 7: Ao = 0.056, 
a = —0.73, and j3 = 1.44. Just as for model C, we assume 
no alignment for L3-5 blue galaxies. We take the 2a lower 
limit for L6 blue galaxies. 

Note that in each model, we used constrained (i.e. a and 
7 fixed) fits to the blue galaxies in each of the luminos- 
ity bins in the r p — (4.7-60)/i _1 Mpc regime. The best- 
fit amplitudes A = wg + (20h~ Mpc) assuming the model 
A scale dependence (a = -0.71) are +0.003 ± 0.026, 
+0.017 ± 0.022, -0.015 ± 0.035, and +0.27 ± 0.25 Mpc 
for the L3, L4, L5, and L6 blue samples respectively, where 
95 per cent confidence errors are given. For model B/C/D 
scale dependence (a — —0.73), we find very similar re- 
sults: w 6+ (20h~ 1 Mpc) = +0.003 + 0.026, +0.017 + 0.022, 
-0.015 + 0.034, and +0.27 + 0.25 ft" 1 Mpc for the same sam- 
ples. For presentation (Table [6} we have re-scaled all values 
to z — 0.3 using the specified choice of 7 and the mean 
redshifts of the samples (see Table [3J. In each case we have 
used the "L3-5 blue" amplitude for the blue galaxies fainter 
than L3, since we have no better constraint for them (the 
SDSS Main sample does not probe an interesting volume at 
fainter luminosities). 

The models are summarized in Table [6j along with es- 



timates of the fractional GI contamination as quantified by 
the ratio C; G/ /Cp G . For most cosmic shear studies, model 
A can be taken as a "2cr" estimate of the GI contamination. 
Model B should be viewed as a "best guess," with the caveat 
that we have no detection of GI in the L3-5 blue samples. 
If the reason for this nondetection is that the GI signal in 
these samples is much less than the current upper limits, 
then model C may be more realistic than B. 

As an example of these models, we show in Figure [TO] 
the GI contamination for a cosmic shear survey that mea- 
sures galaxies down to R = 24 (median redshift 0.6). The 
results in the figure are not necessarily representative of cos- 
mic shear surveys to the specified depth, since the sources 
will not all be weighted equally (and some may be rejected 
from the analysis due to being poorly resolved, or having 
poorly constrained photo-zs with the particular bandpasses 
used in the survey). The fractional contamination (7; /Cp 
is —33, —6.5, —2.1, and —1.5 per cent for models A-D re- 
spectively at I = 500 (these values are given in Table [6]). 
As can be seen from the figure, these fractions are almost 
independent of I because the intrinsic alignment power spec- 
trum P s #i(k) oc k~ 2 ~ a = fc-i-^io.ie ( 95 %CL) has roughly 
the same fc-dependence as the matter power spectrum. This 
was actually measured for the LRGs; setting the blue galax- 
ies (which have no detection) to the same slope as the red 
galaxies was a modeling assumption. The fraction of the 
contamination coming from blue galaxies is 93, 69, 3.3, and 
0.4 per cent for models A-D. It is clear that refining the 
models for blue galaxies should be a high priority for fu- 
ture work since these dominate the uncertainty in the GI 
estimates (though it is not clear whether they dominate the 
GI contamination). Note that the contamination estimates 
presented are on the power spectrum; the fractional contam- 
ination on the amplitude as will be less by a factor of ~ 2 
since the power spectrum is roughly proportional to erf. In 
particular for the survey to R — 24 at I = 500, the reduc- 
tion in shear power spectrum C; 77 is equivalent to changing 
a 8 by Aa 8 = -0.10, -0.02, -0.005, and -0.004 for models 
A-D, respectively. 

Integrated over all scales, the GI contamination we pre- 
dict in the pessimistic model (A) results in a la error for 
a survey measuring all galaxies to R = 24 with combined 
shape and measurement noise 7 rms = 0.3 over a region of sky 
of only 4deg 2 . Thus if model A is correct then GI contami- 
nation may already be important for the current generation 
of lensing surveys. Models B, C, and D predict an error of 
la for sky coverage of 100, 1000, and 2000 deg 2 respectively; 
thus even in the optimistic cases GI contamination will be 
significant and will have to be removed in future surveys. 

We have also investigated the dependence of the con- 
tamination estimates on the magnitude cut, as shown in Fig- 
ure[TT] For this calculation the COMBO-17 luminosity func- 
tion was also used, hence a cutoff in the redshift distribution 
at z = 1.2 was imposed. The amount of contamination gen- 
erally decreases with survey depth: for our "central" model 
it is twice as severe at i? max ~ 22.5 than at i? max ~ 24. 

We also could have constructed a more extreme pes- 
simistic model than model A by taking the 95 per cent con- 
fidence upper limit in each of bins L3, L4, and L5 separately; 
for the above-mentioned toy cosmic shear survey to R — 24, 
this "AA" model leads to —41 per cent contamination at 
I — 500 instead of —33 per cent. Such a model would of 
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Table 6. The four GI models used here to assess contamination. The power law indices are denned by ws+(r p ) oc (1 + z) 1 . The 
amplitudes are normalized to erg = 0.751; for cosmologies with different values of as, the bias of the density tracer scales as b oc a^ 1 and 
the GI contamination scales as Wg+ oc eg. "Fractional contamination" here is defined by C; G/ /C GG . 
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+0.056(L/L ) 1 ' 58 


-0.420 


-0.332 






blue 


-0.71 


+1.47 


+0.028 (L3-5); +0.61 (L6) 






B 


Central 


reel 


-0.73 


-0.56 


+0.059(L/L ) 1 ' 48 


-0.103 


-0.065 






blue 


-0.73 


-0.70 


+0.005 (L3-5); +0.25 (L6) 






C 


Optimistic 


red 


-0.73 


-0.56 


+0.059(L/L ) 1 - 48 


-0.042 


-0.021 






blue 


-0.73 


-0.70 


(L3-5); +0.25 (L6) 






D 


Very optimistic 


red 


-0.73 


-3.29 
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Figure 10. Upper panel: The cosmic shear power spectrum (solid 
line) and GI power spectra (dashed lines) for all galaxies with 
R < 24. We plot the absolute value of GI since it is nega- 
tive for all models. Lower panel: The fractional contamination 
of the power spectrum \Cf II \/Cf l (solid lines), compared to the 
change produced by several changes in the cosmological param- 
eters |ACp G |/C, GG (dashed lines). This ranges from ~ 33 per 
cent for model A to ~ 1.5 per cent for model D. We have plotted 
the range out to I = 2000 (k ~ 1/iMpc" 1 at the typical source 
redshift) since at smaller scales we cannot convert from w g ^.(r p ) 
to wg+(r p ) (see Section [jjl. For the variations of cosmological pa- 
rameters, f2 m is varied at fixed Cl m h 2 and n s is varied at fixed 
o"8- 
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Figure 11. The GI contamination at I = 500 for the four models 
as a function of survey depth. We considered all galaxies with 
i?-band magnitude less than the specified cutoff, and z < 1.2. 



course be unrealistically pessimistic for surveys covering a 
broad range in galaxy luminosity, since by taking the 2a up- 
per limit in each luminosity bin and averaging these numbers 
we obtain a > 2a upper limit on the average. 



9.2 Constraining as 

Finally we come to the issue of greatest practical importance 
in the near term: how should one account for GI contami- 
nation in as measurements from cosmic shear? We have al- 
ready seen that if our more pessimistic models are correct 
then the contamination may be significant compared to the 
uncertainties in some of the recent measurements. In this 
case it is essential to correct for the GI effect and include 
the range of allowed GI models in determining the error 
bars on as ■ Since the present surveys are measuring a single 
amplitude as (with some dependence on other cosmological 
parameters such as O m ) what we really need is a probability 
distribution for Acrg (which also may depend on other pa- 
rameters) over which one can marginalize. In future surveys 
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that measure several parameters we will need a multivari- 
ate distribution. Therefore we view the model here as a first 
step: it leaves much room for future work to reduce the un- 
certainties in Aug, particularly to pin down the several-tr 
tails of the distribution (which are treated rather crudely 
here), and to extend the model to incorporate multiple pa- 
rameters. 

In constructing this distribution for Acts, we first note 
that the error bar on the GI contamination is highly asym- 
metric due to the nature of the redshift extrapolation: if the 
redshift exponent 7 is very negative or zero there is very lit- 
tle contamination, but if it is positive then it is possible to 
have very large contamination. In all cases the central model 
predicts relatively small contamination (a few per cent) but 
the pessimistic model based on the 2a upper limit to 7 could 
be a factor of several worse, and the optimistic models are 
a factor of several better. Therefore we recommend that for 
measurements of as one should marginalize over contamina- 
tion Ao"8 with a lognormal probability distribution. To be 
more explicit, we recommend the following prescription: for 
the distribution of galaxy luminosities, redshifts, and colours 
in a particular cosmic shear survey, compute the pessimistic 
(A) and central (B) power spectra C; GI (A, B). Then compute 
the induced error Aas (A, B) for these two models. One can 
then write 



as (observed) = erg (true) — x, 
where x is lognormally distributed: 



P(x) 



2ir ax 



exp 



1 



hi- 



(25) 



(26) 



for x > and P(x) = for x < 0. The median of this 
distribution is the central model, i.e. xo = — Acts(B), and 
the standard deviation is chosen to place the pessimistic 
model at 2a: 

Usually one uses a Markov chain Monte Carlo (MCMC) 
method to estimate cosmological parameters from a com- 
bination of data sets; in this case one should marginalize 
over x with the prior given by Eq. (|26[1 . This could be done 
by including x as a nuisance parameter in the MCMC, or 
(probably faster) by including the integral over x as part of 
the cosmic shear likelihood function. 

Note that this method has not explicitly used the op- 
timistic models in constructing the probability distribution 
for x, assuming instead that the distribution in In a; is sym- 
metric. In practice this is probably not a serious deficiency 
for two reasons. First, for the cases we have investigated, 
the ratio of contamination for model A to model B is in- 
deed similar to that of model B to model D (5.1 vs. 4.3 for 
R < 24, I = 500; 4.1 vs. 3.0 for R < 23, I = 500), that is 
the error in "log contamination" is close to symmetric. Sec- 
ondly, the error in erg in the optimistic and very optimistic 
models is usually < 0.01, which is negligible compared to 
the purely statistical errors from the current generation of 
weak lensing surveys. This means that while the details of 
the pessimistic tail of the distribution matter (because they 
affect whether cosmic shear can rule out high erg), the de- 
tails of the optimistic tail do not - the optimistic tail might 
as well be "piled up" at zero contamination. 



We believe the method described here is adequate for 
the current generation of cosmic shear surveys in which the 
cosmological constraint is a single amplitude, generally re- 
ported as a value of as with some dependence on Q m . How- 
ever, it will not be adequate for future surveys that will mea- 
sure the redshift and scale dependence of the signal, in which 
the GI contamination must be described by more than one 
number. Also these surveys will push the statistical errors 
on the cosmic shear signal to the < 1 per cent level, i.e. ac- 
cording to the models presented here the GI contamination 
will dominate over statistical uncertainty. In order to make 
use of this data, we will need more external information to 
better constrain the GI models, or use internal information 
from the cosmic shear surveys themselves to simultaneously 
constrain the pure lensing signal, the GI signal, and (if ap- 
plicable) t he II signal. Strategies for this are discussed in the 
litera ture (|Hirata fc Seliakll2004 lKinsll2005l : iBridle fc Kingl 
12007ft . 

In summary, we have constructed a basic model for the 
GI contamination that can be practically integrated into 
current and near- future lensing constraints on as . It is a min- 
imal model, making the simplest assumptions in some cases 
(e.g. power law scale, luminosity, and redshift dependence). 
Nevertheless, it fulfills the basic criteria of being consistent 
with the data and covering the range of allowed values of 
the most important uncertain parameter (7). Aside from 
improving the statistical uncertainties on GI model param- 
eters, future work covering a range of redshifts would also 
improve the robustness of the model by enabling us to con- 
strain more model parameters, such as deviations from the 
power law ~ (1 + z) 1 or different values of a and 7 for blue 
and red galaxies. 



10 CONCLUSIONS 

In this work, we have attempted to further characterize the 
correlation between the intrinsic shear and the density field 
on both small and large scales. To this end, we have char- 
acterized the strength of this alignment as a function of 
galaxy transverse separation, luminosity (focusing on the 
bright end, where it is most prominent), colour, and red- 
shift to z ~ 0.6. In addition, we have established that the GI 
amplitude of LRGs is not significantly different for BG and 
non-BG galaxies for scales above the virial radius, though 
we were unable to verify the degree to which this state- 
ment is true for cases of large luminosity gap. We found 
that the shape-density alignment signal for red galaxies in- 
creases strongly as a function of the luminosity. The scale 
dependence w g +(r p ) oc r -°- 73 ±° 17 ; s very similar to that 
of the matter power spectrum, which is expected if galaxy 
alignment is determined by the local tidal field. Our findings 
that the density-shape alignment effect is present for non- 
BGs up to large scales implies that our original hypothesis 
that the effect is due only to BGs needs to be revised. 

There is no detection of GI correlation for blue galax- 
ies, with the possible exception of the very brightest bin 
(L6 blue). In this case there is a ~ 2.4cr detection if one 
fixes the scaling with r v to the value for the red galaxies. If 
one does a "blind" search over power laws (as done in 
Section f5.ip the signal is not statistically significant. When 
combined with the fact that we have no detection in the 
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lower-luminosity blue bins which have many more galaxies, 
we believe that the L6 blue GI signal is not robust, and ob- 
servationally many of these galaxies are quite close to the 
blue versus red division. In any case the L6 blue galaxies 
are very rare (despite their luminosity, they comprise < 1 
per cent of the flux-limited SDSS Main sample!) and their 
measured GI signal would contaminate the toy cosmic shear 
survey considered in Section [9] at the < 0.1 per cent level 
if taken at face value. Much more important for the cosmic 
shear programme are the numerous fainter (L ~ L*) blue 
galaxies. The GI signal for these objects is consistent with 
zero, however if the signal is near our upper limits, they may 
affect the observed shear power spectrum by of order 10 per 
cent for surveys to R — 24. Reducing the uncertainties in 
the measurements for the L ~ blue galaxies should be a 
priority for future work. It would also be desirable to learn 
more about the GI correlation wg+(r p ) at very small scales 
where linear galaxy biasing (which underlies the method- 
ology of Section [7} breaks down; this would likely involve 
both the w g +(r p ) measurements presented in Appendix lAl 
and simulation or halo model results on the relation between 
galaxies and mass at small scales. 

We have presented several models of the GI signal in 
Section [9] spanning the range from pessimistic assumptions 
(model A) to optimistic (model D). We believe these models 
will be useful for others in comparing against other datasets 
and/or estimating levels of contamination in various surveys. 
In particular it would be useful for lensing surveys to repeat 
the calculations in Section [9] taking into account their ac- 
tual redshift /luminosity /colour distribution instead of the 
toy distributions used here. 

These results should also be useful in comparing against 
simulations to determine the physical cause behind this ef- 
fect, especially in conjunction with previous results on II 
contamination and future work on contamination of three- 
point functions, which may give additional power to dis- 
criminate between intrinsic alignment models. A fuller un- 
derstanding of the physics behind these effects may allow us 
to extend these fitting relations to higher redshift than is 
currently allowed by the data. 

Methods have been proposed (e.g., King 2005) to re- 
move intrinsic alignment contamination using parameterized 
GI and II correlation models. Our determination of the red- 
shift evolution and luminosity scaling of GI correlations will 
help make these methods more feasible in practice. 

Finally, the GI models presented here (and any im- 
proved models that incorporate future observations) may 
be used to determine criteria for excluding galaxies from fu- 
ture cosmic shear surveys to obtain a galaxy sample with 
the lowest possible level of intrinsic alignment contamina- 
tion. Surveys with imaging data in multiple bands should 
fairly easily be able to remove the bright red galaxies that 
seem to show the strongest GI contamination. Future work 
with simulations will be necessary to determine the efficacy 
of this plan. This plan is most likely to be effective if the blue 
galaxy contamination is found in future work to be negligi- 
bly small; if this is not the case, then m ore sophisti cated 
methods to remove GI such as templates ([Kind 20051) and 
separ ation based on redshift dependence l|Hirata fc Seliakl 
120041 ') will be essential to realising the promise of the cosmic 
shear programmes. 
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APPENDIX A: GALAXY DENSITY-SHAPE 
CORRELATION FUNCTION DATA 

This appendix lists the correlation functions w g +(r p ) for the 
samples used in our fits, and their correlation coefficients. 
The table headers describe the subsamples used to trace 
the intrinsic shear field; the galaxy density ( u g" ) was traced 
using the full sample (i.e. SDSS Main, SDSS LRG low-z, 
SDSS LRG high-z, or 2SLAQ). It is the bias of the latter 
(described in Section that should be used to convert to 
u>s+(r p ). Table IXT1 lists the measured correlation functions 
and their error bars, while Table IA2I shows the correlation 
coefficients p[ws+(r p ), ws+{r' p )] between different radial bins 
for the same sample. The innermost bin from the 2SLAQ- 
bright sample did not have enough data to determine a boot- 
strap error, so it and its correlation matrix elements are 
marked with a "*" in these tables. 



APPENDIX B: RELATION OF LENSING 
SHEAR-INTRINSIC SHEAR AND 
DENSITY-INTRINSIC SHEAR CORRELATION 
FUNCTIONS 

The purpose of this appendix is to relate the lensing shear- 
intrinsic shear correlation function (7 G (z 3 ) ■ 'y 1 (z\ ))e mea- 
sured in simulations by iHevmans et al.l l|2006bl ) to the 
density-intrinsic shear correlation function ws+(r p ) consid- 
ered in this paper. In particular we would like to know 
whether the lHevmans et all (|2006bh are consistent with our 
observations. In this appendix, we will use z 3 for the source- 
plane redshift and z\ for the lens-plane redshift (i.e. the red- 
shift at which GI contaminati o n is be ing assessed) for con- 
sistency with IHevmans et al.l l|2006bl ); how ever for consis- 
tency w ith our paper we will write -y 1 where IHevmans et al.l 
(2006b]) would write e [defined by their Eq. (4); note that 
the ellipticity components before correction for the shear 
responsivity, which w e denote e*, are denoted by in 
IHevmans etafl (|2006bl )]. 

Our conversion proceeds in two steps: first from 
(■y G (z s ) ■ 7 7 (jzi))e to (n(z s )'y I + (zi))g, and then from 



(ft(z s )7+(0i))e to ws+(r p ). The first step is a straightfor- 
ward but tedious calculus exercise, since 7° and k are dif- 
ferent derivatives of the lensing potential and hence are sim- 
ply related in Fourier space; the second step is trivial and 
involves the usual convergence to surface density conversion 
from galaxy-galaxy lensing studies. 

For our first step, we recall that the shear components 
+ , x can be easily converted to E and B-modes in Fourier 
space by 

7?(«) ± i 7 ?(0 = e 2m) h%(l) + i7l(01. (Bl) 

where cj>(l) — a,rcta,n(l y /l x ) is the position angle of I, and 
a similar formula is written for the intrinsic shear. The 
convergence-density correlation function is then the Fourier 
transform of the cross-power spectrum, 

<k(0)7|(o)> = / -J^rvv^) 

f°° Idl 1 
= -J o MW). (B2) 

Here we have assumed parity invariance, i.e. that the con- 
vergence is correlated only with the S-mode of the intrin- 
sic shear. (The latter may have a B-mode but this cannot 
be correlated with a scalar.) We have also chosen to he 
along the x-axis (i.e. the axis along which stretching cor- 
responds to positive +- component shear) and used the in- 
tegral representation of the Bessel function (Eq. 8.411 of 
iGradshtevn fe Rvzhiklll99i ). A similar result holds for the 
7 G 7 7 correlation: 

/°° / d/ / 
— cr Mm; (B3) 

once again we have used only the £-mode since the lens- 
ing shear 7° is derived from the scalar density field and 
has only _E-modes (aside from small corrections due to 
non-weak shear or multiple deflections). We have replaced 
C T T — * Of since 7°. and k have equal numerical val- 
ues mode by mode (they are both second derivatives of the 
lensing potential), and we have a Jo function instead of — J2 
because there is no factor of e 2 "^^ in this equation [the dot 
product of the lensing and intrinsic shears is unaffected by 
the angle-0(Z) rotation of Eq. (|B1[) ]. 

The conversion between Eq. (|B2[I and Eq. ()B3|I proceeds 
as follows. The cross-power spectrum C^~* can be obtained 
from Eq. (IB3|l by the Hankel transform pair, 

1 r°° 

Cf 7 = 2n (7 G (tf) ■ 7 / (0)}J (M)tfdtf, (B4) 
Jo 

from which we find 

(k(0)7+(O)> = - / <7 G W-7 / (O))e(0,^d 1 9, (B5) 
Jo 

where the Green's function is 

poo 

g{0,&)= dllJo(W)J 2 (lO). (B6) 
Jo 

This Green's function has a closed-form analytical expres- 
sion that can be obtained using the recursion, differentiation, 
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Table Al. The correlation functions w g +(r p ) for each sample. The first column shows the range of r p , and the remaining columns show 
w g +(r p ) and its la uncertainty (i.e. square root of diagonal covariance matrix element). Units are h~ 1 Mpc for all columns. 





max 


Main. L3. blue 


Main. L4. blue 


Main. L5. blue 


Main. L6. blue 


Main. L3. red 


Main.L4.red 


0.3 


-0.6 


0.145 ± 0.198 


-0.177 ±0.194 


0.494 ± 0.443 


-1.127±6.011 


0.117 ±0.330 


0.032 ±0.212 


0.6 


-1.2 


-0.081 ± 0.129 


0.118 ±0.108 


-0.106 ±0.215 


1.604 ± 1.854 


0.259 ±0.147 


0.070 ±0.100 


1.2- 


-2.4 


-0.038 ± 0.069 


-0.020 ± 0.067 


0.100 ±0.134 


1.588 ±0.761 


0.003 ±0.101 


-0.005 ± 0.053 


2.4- 


-4.7 


0.033 ± 0.035 


-0.012 ±0.035 


0.018 ±0.088 


0.590 ± 0.653 


0.053 ±0.057 


0.008 ± 0.053 


4.7- 


-7.5 


0.017 ±0.034 


0.061 ±0.036 


-0.030 ± 0.060 


0.948 ± 0.395 


0.103 ±0.048 


0.024 ± 0.036 


7.5- 


11.9 


0.059 ± 0.042 


0.016 ±0.030 


-0.040 ± 0.050 


0.585 ± 0.329 


0.110 ±0.042 


0.036 ± 0.026 


11.9- 


-18.9 


-0.000 ± 0.032 


0.036 ± 0.026 


-0.048 ± 0.046 


0.227 ±0.313 


0.073 ± 0.041 


0.014 ±0.024 


18.9- 


-29.9 


-0.010 ± 0.025 


0.015 ±0.021 


-0.000 ± 0.045 


0.177 ±0.214 


0.059 ± 0.040 


0.011 ±0.018 


29.9- 


-47.4 


0.003 ± 0.020 


0.007 ±0.021 


-0.007 ±0.029 


0.176 ±0.155 


0.010 ±0.034 


0.008 ± 0.013 


47.4- 


-59.6 


0.017 ±0.026 


0.003 ±0.019 


0.007 ± 0.034 


-0.046 ±0.164 


-0.044 ± 0.035 


0.010 ±0.019 


fj) ? rn i ii _ 


"^p,max 


Main. L5. red 


Main. L6. red 


LRGl,low-2 


LRG2,low-z 


LRG3,low-2 


LRGl,high-z 


0.3- 


-0.6 


0.744 ±0.224 


7.676 ±1.387 


5.075 ±9.556 


1.365 ±8.740 


2.975 ± 6.866 


-1.277 ±8.148 


0.(5 


-1.2 


0.278 ± 0.127 


2.906 ± 0.644 


3.590 ± 1.753 


1.007 ± 2.188 


1.714 ±3.181 


2.784 ± 1.583 


1.2 


-2.4 


0.212 ± 0.087 


0.956 ±0.288 


1.978 ±0.873 


-0.216 ±0.606 


1.633 ±0.769 


0.370 ± 0.526 


2.4 


-4.7 


0.149 ± 0.051 


0.619 ±0.158 


1.033 ±0.327 


0.105 ±0.312 


0.957 ±0.302 


0.693 ± 0.242 


4.7- 


-7.5 


0.039 ± 0.036 


0.420 ±0.132 


0.270 ±0.189 


0.182 ±0.187 


0.362 ±0.190 


0.441 ±0.191 


7.5- 


11.9 


0.062 ± 0.028 


0.243 ±0.107 


0.274 ±0.130 


0.078 ±0.164 


0.708 ±0.135 


0.388 ±0.133 


11.9- 


-18.9 


0.046 ± 0.030 


0.306 ± 0.076 


0.339 ±0.100 


-0.001 ±0.100 


0.250 ± 0.087 


0.349 ± 0.076 


18.9- 


-29.9 


0.092 ± 0.027 


0.239 ± 0.063 


0.135 ±0.063 


0.035 ±0.059 


0.209 ± 0.062 


0.209 ± 0.071 


29.9- 


-47.4 


0.036 ± 0.018 


0.120 ±0.060 


0.053 ±0.050 


0.106 ±0.055 


0.164 ±0.043 


0.138 ±0.035 


47.4- 


-59.6 


0.010 ±0.016 


0.016 ±0.038 


-0.036 ± 0.057 


0.115 ±0.045 


0.092 ±0.041 


0.107 ±0.047 


^p,miii" 


"?"p,max 


LRG2,high-z 


LRG3,high-z 


2SLAQ-faint 


2SLAQ-bright 






0.3 


-0.6 


-7.429 ± 7.872 


21.379 ±9.708 


-3.713 ±8.009 


* 






0.(5 


-1.2 


4.603 ± 3.609 


7.346 ±3.589 


-3.235 ±3.722 


5.380 ±3.223 






1.2 


-2.4 


2.940 ± 1.125 


2.634 ±0.792 


0.425 ± 1.327 


2.717 ± 1.287 






2.4 


-4.7 


1.950 ± 0.444 


1.150 ±0.406 


-0.155 ±0.499 


0.598 ± 0.621 






4.7- 


-7.5 


0.629 ± 0.296 


0.890 ± 0.308 


0.506 ± 0.385 


-0.230 ± 0.396 






7.5- 


11.9 


0.518 ± 0.208 


0.539 ±0.155 


-0.335 ±0.333 


0.469 ± 0.214 






11.9- 


-18.9 


0.592 ± 0.113 


0.365 ±0.107 


-0.168 ±0.161 


0.087 ±0.155 






18.9- 


-29.9 


0.344 ±0.103 


0.290 ± 0.084 


0.004 ±0.172 


-0.020 ±0.138 






29.9- 


-47.4 


0.094 ± 0.080 


0.257 ±0.054 


-0.085 ±0.157 


0.209 ±0.116 






47.4- 


-59.6 


0.103 ±0.077 


0.144 ±0.060 


-0.129 ±0.207 


0.112 ±0.143 







and orthonormality relations for Bessel functions: 



-Jo(W) + -Ji(W) 



&UJo(W) 

dllJ (W)J ^>) -!- j t j >\l .l,Ah>).J t ih)) 
2 • Al 

dw 



5(0- 


0) 


9 

5(0- 


*) 


9 

5(0- 


0) 



Aw Ui(lw)Ji(l0) 
5(w - 0) 



+ 3a e (*-*)' 



(B7) 



where O is the Heaviside step function. (The third line has 
used the fact that — Ji is the derivative of Jo.) Therefore 



<*(0)7+(O)> 



-( 7 G (e). 7 J (0)) 



(P 



(7 G W ■ 7 7 (0))i?dtf. 



(B8) 



iHevmans efall (|2006bl ) fit their results to the functional 
form 

A£ 



(~ 1 G ( Zs )- 1 I (z 1 ))e = 



(B9) 



where A and 9o are free parameters and £ = DiDi s /D 3 is 
the lensing strength. Plugging this into Eq. (|B8|I yields 

Our second step is to convert this into a density-intrinsic 
shear correlation function. To do this, we recall the Born 
(single-deflection) approximation for the convergence, 



(Bll) 



where 5 is the density perturbation. If the correlations be- 
tween k(z s ) and the intrinsic shear field at a particular red- 
shift zi arise near zi , then we may change variables from dis- 
tance from the observer r to radial separation II = r— r(zi), 
multiply both sides by the intrinsic shear field, and take the 
average: 

^(z a W + ( Zl )) e = j dn£ a+ (r PJ II), (B12) 

where r p — rO is the separation and we can pull £ (which 
varies slowly with lens redshift) out of the integral because 
only a small range at z ~ zi contributes. The last integral is 



© 0000 RAS, MNRAS 000, 000-000 



24 Hirata et al. 



Table A2. Correlation matrices for the data shown in Table IA1I The units are per cents (i.e. all correlation matrix elements have 
been multiplied by 100). In each of the 4 blocks of this table, there is an upper (above-diagonal) triangle, which applies to the sample 
designated "upper," and a lower (below-diagonal) triangle, which applies to the sample designated "lower." Note that the diagonal 
correlation coefficients are equal to unity (100 per cent). 
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ws+(r p ), so combining with Eq. f|B10f) we have 
t ^ 2(?A ( 1 2 20o, + e o \ 

Ws+ ^ ) = 3^{eTe-o-e + -^ ln -er)- (B13) 

This is the equation we have used to compare the 
iHevmans et al.l (|2006bl ) results to ours in Section [8.21 
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